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INTRODUCTION 

Calculations  have  been  reported  from  the  Naval  Postgraduate  School  for  the  results  of 
adiabatic  internal  explosions  in  systems  containing  a  variety  of  C-H-N-0  fuels  in  air  with  or 
without  the  addition  of  the  active  metals  magnesium  or  aluminum.1 **7  The  present  study 
reflects  the  extension  of  the  studies  to  titanium  as  the  active  metal. 

As  before,  metal  and  fuel  are  considered  introduced  into  a  constant  volume  of  1  m3 * 5 *, 
initially  at  25°C.  Products  are  assumed  distributed  evenly  throughout  the  volume.  The 
process  is  treated  as  adiabatic,  and  the  ideal-gas  approximation  is  assumed  to  hold. 

The  fuels  considered  are  listed  in  Table  1,  along  with  a  few  pertinent  properties. 
Oxygen  balance  is  computed,  on  a  mass-percent  basis,  as  the  excess  (deficiency,  when 
negative)  of  oxygen  in  the  fuel  relative  to  the  production  of  water,  carbon  monoxide,  and 
elemental  nitrogen. 

In  all  the  following,  the  symbol  C  will  represent  the  charge,  or  mass  of  fuel,  per  cubic 
meter  of  air  and  M,  the  mass  of  metal  (titanium)  per  cubic  meter.  Charge-to-metal  {C/M) 
ratios  and  total  concentrations  (C  +  M,  in  kilograms  per  cubic  meter)  were  varied  over  the 
range  0.1  to  10  for  each. 

In  several  of  the  tables,  a  code  is  used  to  indicate  the  condensed  phases  present:  L  = 
liquid  oxide  solution;  C  =  TiC;  N  =  TiN;  G  =  graphite;  A  =  Ti30&;  B  =  T1O2. 


1  Naval  Weapons  Center.  Peak  Overpressures  for  Internal  Bluet ,  by  G.  F.  Kinney,  R.  G.  S.  Sewell,  and  K.  J. 
Graham.  China  Lake, Calif.,  NWC,  June  1979.  (NWC TP 6087,  publication  UNCLASSIFIED.) 

^  Naval  Weapons  Center.  Reactive  Metals  in  Internal  Explosions:  The  Combustion  of  Magnesium  in  Air ,  by 
R.  A.  Reinhardt.  China  Lake.  Calif.,  NWC,  February  1978.  (NWC  TM  3429,  publication  UNCLASSIFIED.) 

3  Naval  Weapons  Center.  Adiabatic  Computation  of  Internal  Blast  for  Magnesium -Cased  Charges  in  Air,  by 
R.  A.  Reinhardt.  China  Lake,  Calif.,  NWC,  April  1979.  (NWC  TM  3820,  GIDEP  E202-I481,  publication 
UNCLASSIFIED.) 

*  Naval  Weapons  Center.  Adiabatic  Computation  of  Internal  Blast  for  Aluminum -Cased  Charges  in  Air,  by 
R.  A.  Reinhatdt  and  A.  K.  McDonald.  China  Lake,  Calif.,  NWC,  January  1982.  (NWC  TP  6287,  publication 
UNCLASSIFIED.) 

5  Naval  Weapons  Center.  A  Working  Model  for  the  System  Alumina-Magnesia,  by  R.  A.  Reinhardt.  China 
Lake, Calif.,  NWC,  May  1983.  (NWC  TP6433.GIDEP  E392-0754,  publication  UNCLASSIFIED.) 

3  Naval  Weapons  Center.  Computer  Program  for  Internal  Aluminum-Fuel- Air  Explosions,  by  R.  A. 
Reinhardt.  China  Lake.Calif.,  NWC.  May  1983.  (NWC  TP6449,  GIDEP  E413-0319,  publication  UNCLASSIFIED.) 

7  Naval  Postgraduate  School.  Internal  Exploeions  of  Reactive  Aluminum  with  a  PBX  in  Air,  by  R.  A. 
Reinhardt.  Monterey,  Calif.,  NPS,  August  1983.  (NPS-61-83-011-PR,  publication  UNCLASSIFIED.) 
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TABLE  1.  Properties  of  the  Fuels. 


Designation 

Name  and  chemical  formula 

AUf,  298. 
kJ/mole 

Oxygen  balance, 
to  CO  +  H20,% 

PETN 

i 

Pentaerythritoltetranitrate, 

C5H8N4Oi2 

-489.8 

+  15 

NC 

Nitrocellulose,  13.3%  N; 
C6h7N2.5°10 

-812.9 

+  3 

HMX 

sym-Cydotetramethylene- 
tetranitramine,  C4HgNg08 

+  77.4 

0 

Pentolite 

50%  PETN,  50%  TNT; 

C6.I6H6  25N3.4I08.5 

-97  9 

-5 

Comp  8 

65%  RDX,  35%  TNT; 

Cl  96H2.53N2.22°2  68 

+  114? 

-9 

TNT 

2,4,6-Trinitrotoluene, 

C7H5N3O6 

+  30.7 

-25 

n2h4 

Hydrazine 

-165.0 

-100 

c2h4o 

Ethylene  oxide 

-47.29 

-109 

Carbon 

Graphite 

0 

-133 

C6Hi4 

Hexane 

+  100.31 

-242 

Internal  energies  of  formation  are  per  mole  of  formula  indicated. 
Data  were  taken  or  computed  from  Ref.  1 


BASIS  OF  CALCULATIONS 

As  is  pointed  out  in  Ref.  4,  the  adiabatic  restriction  requires  that  a  temperature  be 
found  such  that  the  sums  of  the  internal  energies  of  the  equilibrium  mixture  of  products  at 
that  temperature  must  equal  the  initial  internal  energy  of  formation  of  the  chosen  fuel.  The 
pressure  is  then  computed  from  the  ideal  gas  law,  taking  into  account  the  total  number  of 
moles  of  gas  present  at  equilibrium  at  the  temperature  The  overpressure  is  found  by 
subtracting  1  bar. 

As  before,4  thermochemical  data  are  represented  by  the  five-parameter  expression 
U  -Bi  +  B2T  +  B3T2  +  B4lnT  +  BSIT 

where  T  is  the  Kelvin  temperature.  In  Table  2  are  listed  the  internal  energy  parameters  for 
the  25  gaseous  and  15  condensed-phase  species  considered.  Those  for  titanium  species  were 
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TABLE  2.  Internal  Energies  of  Combustion  Products. 
Expressed  as  a  function  of  x  =  T/1000  (in  kilokelvins); 
values  given  in  joules/mole:  U(T)  =  Bi  +  B2x  +  B3T2  +  B^nx  +  B5/X 


Substance 

Bi 

82 

83 

84 

85 

Ti 

310931 

11660 

75 

1364 

664 

TiO 

2640 IS 

30151 

272 

-2961 

1040 

TiQ2 

-  55866 

60797 

-355 

-22628 

-4116 

Ti  + 

1029558 

29682 

-252 

-49343 

Ar 

-3718 

12473 

0 

0 

CO 

-84115 

31136 

-25 

-6218 

C02 

-361013 

56789 

67 

-9201 

H 

213029 

12473 

0 

0 

OH 

166505 

34786 

-88 

-20093 

”2 

136535 

32890 

313 

-21058 

h2o 

50545 

63224 

-395 

-  45949 

NO 

105501 

30619 

10 

-4353 

n2 

38785 

31651 

-61 

-7831 

0 

214445 

9604 

298 

4428 

°2 

32105 

29476 

662 

-6265 

TiC(l) 

-617142 

62760 

0 

0 

TIC  (S) 

-762144 

27917 

6225 

14807 

Ti  <l> 

-465953 

35564 

0 

0 

T.(s) 

-695970 

-10976 

8134 

35009 

C« 

-40329 

22904 

320 

10 

CN" 

102390 

31551 

-69 

29 

CN 

691554 

38828 

740 

-37192 

C2H 

758476 

66613 

239 

-  43882 

c2n 

538335 

54371 

-15 

-1052 

HCN 

284944 

62760 

-302 

-  26216 

HNCO 

36816 

83592 

-474 

-28281 

HCO 

90276 

56166 

-353 

-  19233 

ch2o 

54534 

85563 

-609 

-32362 

c2h2 

474890 

92019 

-27 

-42593 

C3 

906056 

48907 

282 

-16185 

TlO(l) 

-1040198 

66944 

0 

0 

TiO  (s) 

-1011246 

56480 

4163 

0 

T12O3  <l) 

-  1942190 

156900 

0 

0 

Ti  203(s> 

-2041671 

145104 

2725 

-30 

T13O5  CD 

-2875209 

234304 

0 

0 

Ti305(s) 

-2963674 

174570 

16844 

266 

Ti02  (1) 

-  1395260 

87864 

0 

0 

Ti02  (s) 

-1436271 

62766 

5687 

159 

TiN  (1) 

-  770564 

62760 

0 

0 

TiN  (S) 

-  847427 

38404 

5099 

4361 

Substances  are  gaseous  unless  otherwise  specified. 
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calculated  from  data  in  the  JANAF  tables, 8  recomputing  the  data  so  as  to  accommodate  the 
universal  choice  of  titanium  vapor  as  the  reference  state.  Parameters  for  the  other  species  are 
taken  from  Ref.  4. 


EQUILIBRIUM  CONSIDERATIONS 
THE  SYSTEM  TITANIUM-OXYGEN 

Unlike  magnesium  and  aluminum,  the  metals  studied  earlier,  titanium  forms  a 
number  of  solid  oxides.  The  compounds  TiO,  Ti203  and  Ti02  arc  well  known,  and  their 
formation  accords  with  the  solution  chemistry  of  this  transition  metal.  The  phase  diagram  of 
the  system  Ti-0  is  reasonably  well  established  up  to  the  melting  points.1®  Maxima  in 
melting  point  appear  to  correspond  to  T^O,  Ti^Ai,  Ti3Gj  and  Ti02-  TiO  appears  to  melt 
incongruently.  There  is  extensive  solid-solution  formation  and  considerable  temperature- 
dependent  polymorphic  transition.  Oxygen  is  soluble  in  titanium  metal  to  the  extent  of  about 
25  atom-percent,  yielding  solid  solutions  which  show  metallic  conductance.  In  general,  the 
lower  oxides  are  all  nonstoichiometric  and  show  semiconductor  properties. 

The  phase  diagram  gives  no  direct  evidence  regarding  the  liquid  phase,  but  the 
appearance  is  not  inconsistent  with  that  of  a  system  with  no  discontinuities  in  the  liquid 
phase. 


For  these  reasons,  it  has  been  assumed  in  the  present  study  that  the  liquid  phase  in 
the  system  Ti-0  is  continuous;  that  is,  that  there  is  complete  miscibility  from  the  metal  to  the 
composition  Ti02  This  liquid  is  assumed  to  be  an  ideal  solution  of  the  liquids  Ti,  TiO,  T12O3, 
Ti3Os  and  Ti02-  T^O  was  not  included  since  no  thermochemical  or  equilibrium  data  appear 
to  be  available  for  it. 

The  complexities  of  the  solid  phases  have  been  ignored  in  the  present  study.  Only  a 
small  number  of  points  were  at  temperatures  below  the  melting-point  curve;  most  of  these 
were  for  oxygen-rich  systems  in  which  solid  T1O2  formed  or  for  carbon-rich  systems  in  which 
TiC  but  no  oxide  at  all  resulted.  Pure  solid  phases  were  assumed  in  those  few  cases  where 
solid  Ti305,  alone  or  with  Ti02,  appeared. 

.In  the  computer  program  the  presence  of  "liquid  oxide”  (that  is,  the  liquid  solution  of 
oxygen  in  titanium)  required  in  part  that  the  temperature  be  above  2143  K  (the  melting  point 
of  Ti02  and  the  highest  temperature  on  the  melting-point  curve).  This  criterion  was  adopted 
for  the  sake  of  simplicity,  and  proved  satisfactory  in  all  but  a  handful  of  cases  where  a  liquid 
was  obtained,  but  below  2143  K.  These  points  were  recomputed,  now  waiving  the 
temperature  requirement  for  liquid  oxide.  It  was  found  that  this  change  made  very  little 
difference  (at  most  a  few  tenths  of  a  bar)  in  the  computed  pressure. 


®  National  Bureau  of  Standards.  JANAF  Thermuchemical  Tables ,  2nd  edition,  by  D.  K  Stull  and  H. 
Prophet.  Washington,  D.C.,  NBS,  June  1971 . 

®  R.C.  DeVries  and  R.  Roy,  Am.  Ceramic  Soc.  Bull.,  Vol.  3309541,  pp,  370-72. 

P.G.  Wahlbeckand  P,  W.  GiUea,,/.  Am.  Ceramic  Soc.,  Vol.  49  ( 1966),  pp.  180-83. 
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A  set  of  calculations  was  carried  out  for  the  system  TNT-Ti-air  using  the  quite 
different  (and  less  reasonable)  assumption  that  the  liquid  oxide  phases  were  pompletely 
immiscible  with  one  another  and  with  the  liquid  metal.  Results  were  appreciably  different,  as 
to  both  adiabatic  temperatures  and  pressures  found.  The  pr  essures  tended  to  be  higher  (by  up 
to  12%)  even  though  the  temperatures  were  lower.  This  result  was  in  accord  with  the 
reduction  in  fugacities  accompanying  solution  formation. 

It  should  be  remarked  that  it  was  assumed  that  TiC,  TiN  and  graphite,  when  present, 
existed  as  pure  phases. 


EQUILIBRIUM  DATA 

As  in  the  previous  study,4  equilibrium  constants  of  formation  were  fitted  to  the  four- 
parameter  equation: 

log  10  K  =  At  +  A  2/(A3  +  T)  +  A4T 

In  Table  3  are  listed  the  equilibrium  constant  parameters  for  all  the  species  considered 
(except  for  the  elements  in  the  reference  states,  for  which  all  the  parameters  are  zero).  Those 
for  the  titanium  species  were  computed  from  the  JANAF  tables,8  with  the  data  below  3591  K 
(the  normal  boiling  point  for  titanium)  recomputed  to  take  into  account  the  choice  of  titanium 
vapor  as  the  reference  state,  (This  choice  was  made  to  avoid  the  discontinuities  that  would 
otherwise  occur  at  the  melting  and  boiling  temperatures  of  titanium.)  Parameters  for  the 
other  species  are  taken  from  Ref  4. 

In  application,  each  Kp  (expressed  in  partial  pressures)  is  first  converted  to  a  Kn  (in 
terms  of  mole  numbers).  Then  for  any  species  the  number  of  moles  is  given  as  a  known 
function  of  Kn  and  the  "master  variables":  X  =  VO2;  Y  =  VH2;  Z  =  n/NV.  TL;  and  Acc,  which 
is  the  activity  of  carbon  (standard  state,  graphite).  (Formulas  in  the  last  sentence  refer  to 
number  of  moles.) 


EQUILIBRIUM  CALCULATIONS 

For  argon,  the  mole  number  is  always  0.4036,  the  number  of  moles  in  1  m3  of  air  at 
298  K,  1  bar.  For  each  of  the  Five  remaining  elements,  a  material  balance  equation  may  be 
written  representing  conservation  of  the  number  of  moles  of  atoms  of  each  element.  The  mole 
number  of  each  of  the  chemical  species  present  is  a  function  of  those  of  the  elements  in 
standard  states  and  the  activity  of  carbon;  hence,  five  simultaneous  equations  in  five 
unknowns  (the  five  "master  variables"  referred  to  earlier)  are  obtained. 

_  In  actuality  the  maximum  number  of  unknowns  to  be  considered  is  four,  since  Y  = 
vh2  can  always  be  found  in  closed  form  as  a  function  of  the  others.  Considering  also  the 
phase  rule  restriction,  the  number  of  variables  to  be  solved  for  is  four  less  the  number  of 
condensed  phases.  The  Newton-Raphson  method,  as  described  in  Ref.  4,  is  then  used  to  solve 
this  set  of  simultaneous  nonlinear  equations. 
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TABLE  3.  Equilibrium  Constants  of  Formation  of  Products. 

(Base  1 0  logarithm  of  the  formation  constant  of  the 
indicated  substance,  expressed  as  a  function  of  t  =  T/tOOO 
(in  kilokelvins):  logioK  =  A,  +  A2A  A3  +  1)  +  A4t) 


Substance 

Al 

A2 

A3 

A4 

TiO 

-2.256 

23  248 

-0  018 

-0.133 

T1O2 

-7.510 

37  603 

-0.008 

-0.067 

Ti  + 

4.082 

-39.359 

0.079 

0  027 

CO 

4.593 

6.108 

0.030 

-0.067 

C02 

0.087 

20.642 

0  001 

-0.024 

H 

3.132 

-12.016 

0.019 

0015 

OH 

0.769 

-1  969 

-0.016 

-0015 

H20 

-3.056 

13.305 

0.014 

-0  005 

NO 

0.710 

-4  300 

0.008 

-0  009 

0 

3.497 

-13  439 

0.010 

0  006 

TlC(l) 

-6.352 

29.181 

-0.002 

0.019 

TiC{s) 

-8.124 

33.929 

-0  003 

0  118 

Ti  (1) 

-6  033 

21  383 

-0  060 

-0.006 

Ti  (s) 

-7  550 

24  36! 

0000 

0.165 

CN- 

0826 

0  005 

-2.418 

0  047 

CN 

5.092 

-22  218 

-0012 

-0.016 

C2h 

6.529 

-24  269 

-0.012 

-0.043 

c2n 

7.060 

-29  343 

-OOOO 

-0.066 

HCN 

1.607 

-6  807 

-0,012 

-0  008 

HNCO 

-25.740 

96.491 

2.563 

2  862 

HCO 

2  221 

1  209 

0.219 

-0.060 

ch2o 

-2.036 

6  770 

0018 

-0  016 

CjHj 

2  655 

-11  332 

-0016 

0  000 

C3 

10  850 

-41.395 

-0016 

-0  186 

TiO  (1) 

-8.512 

44816 

-0  050 

0  052 

TiO  (s) 

11.560 

50  682 

-0.009 

0  348 

T12O3  (1) 

-22.269 

112  701 

-0.054 

0  030 

-28.610 

126.047 

-0  007 

0  607 

-36.762 

182  554 

-0.039 

0  085 

S11E2  SM 

-45  523 

200  007 

-0.003 

0  984 

tio2  (i) 

-13.912 

66  399 

-0.042 

0042 

T1O2  (s) 

-17.167 

73  460 

-0  003 

0  290 

TiN  (1) 

-5  271 

20076 

-0.672 

-0.437 

TiN  (s) 

-12.371 

41  581 

-0  007 

0  196 

Substances  are  gaseous  unless  otherwise  specified. 
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An  initial  approximation  scheme  is  available  wherein  it  is  assumed  that  oxygen  is 
taken  up  in  the  order  CO,  TiOx,  H20,  COj,  With  insufficient  oxygen  to  proceed  beyond  TiO, 
TiN  is  assumed.  With  excess  carbon,  TiC  or  C  (solid)  is  assumed.  Based  on  these 
assumptions,  initial  values  of  all  the  master  variables  can  be  found  and  the  computation  may 
be  started.  Alter  the  first  computation  in  a  series  it  is  a  matter  of  operator  choice  whether  to 
repeat  the  initial  approximation  or  to  use  the  values  of  the  master  variables  from  the  last  run. 

In  any  case,  after  convergence  of  the  Newton's  method  of  calculation,  tests  are 
performed  for  the  presence  or  absence  of  each  condensed  phase.  If  a  change  from  those 
assumed  has  occurred  the  computation  is  rerun.  This  procedure  continues  until  all  the  tests 
are  satisfied  and  no  change  in  condensed  phases  is  predicted.  Further  detail  on  the 
computational  method  with  a  copy  of  the  computer  program  is  given  in  Appendixes  A  and  B  to 
this  report. 


RESULTS 

For  all  the  systems  studied  the  following  were  found:  adiabatic  temperature, 
overpressure,  and  yield  of  each  product  in  terms  of  moles  and  also  as  partial  pressures  for  the 
gaseous  products.  Total  concentration  (metal  plus  fuel)  was  varied  regularly  over  the  range 
0.1  to  10  kg/m3  and  the  fuel-to-metal  ratio  likewise  was  varied  over  the  range  0.1  to  10. 

In  Table  4  the  following  data  are  given  for  titanium  plus  air  (in  the  absence  of  fuel): 
overpressures  (bars),  adiabatic  temperatures  (K),  and  product  yields,  expressed  as  mole- 
percent  for  the  gases  and  total  mole  number  for  condensed  phases.  The  composition  of  the 
liquid  oxide  phase  is  given  as  x  in  the  empirical  formula  TiOx  (x  ranging  from  0  to  2). 


TABLE  4.  Combustion  of  Pure  Titanium  in  Air: 
Overpressure.  Adiabatic  Temperature,  and  Product  Yield. 


Property/ 

system 

Concentration,  kg/tn3 

0.1 

0  2 

04 

1  0 

20 

40 

100 

Overpressure,  bars 

7  3 

9  7 

117 

14.3 

18.4 

26.3 

49  8 

Temperature.  < 

2602 

3423 

4024 

4338 

4393 

4542 

4869 

Mole-% : 

Ar 

1  05 

1  08 

1  07 

0.95 

0.76 

0  56 

0.32 

Ti 

0  01 

7.22 

30  08 

50  59 

72.46 

TiO 

0  12 

486 

17.42 

9  85 

5.22 

2.15 

TiOj 

0  01 

008 

NO 

2.47 

5  05 

409 

005 

001 

N2 

80  84 

81  46 

81  12 

74  26 

59  59 

43.63 

25.06 

O 

0  32 

3  64 

6.54 

0.10 

0  01 

0.01 

15.31 

8  64 

2  23 

...  | 

Liquid  oxide: 

Moles 

2  04 

2  77 

4  55 

10  40 

20  55 

43  2 

1150 

x  in  T-0* 

1  99 

1  79 

1  40 

0  91 

0.57 

0.31 

0  12 

Jib 
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Table  5  gives  overpressures  in  bars,  with  the  temperature  in  parentheses,  for  the 
titanium  air-fuel  mixtures.  Each  portion  of  the  table  is  devoted  to  a  particular  fuel,  and  the 
fuels  are  arranged  in  decreasing  order  of  oxygen  balance,  starting  with  the  most  oxygen-rich, 
PETN,  and  ending  with  the  most  oxygen-deficient,  hexane.  In  each  case,  C  +  M  represents 
the  total  concentration  in  kilogram/meter3,  and  C/M  represents  the  mass  ratio  of  fuel  to 
metal.  For  each  table  entry  a  code  indicates  the  condensed  phases  present;  this  code, 
described  earlier,  is  also  explained  as  a  footnote  to  the  table. 


TABLE  5.  Overpressure  (in  bars).  Adiabatic  Temperature  (in  K)  (in  parentheses), 
and  Condensed  Phases  for  the  Combustion  of  Titanium  With  Fuels  in  Air. 


a.  Fuel;  PETN. 


NWCTP6544 


TABLE  5.  (Contd.). 


r/M 

P.. 

0.1 

0.2 

0.4 

1.0 

2  0 

4.0 

10.0 

b. 

Fuel:  Nitrocellulose. 

7.1 

(2S08) 

L 

6.9 

(2427) 

L 

6.6 

(2294) 

L 

6.1 

(2089) 

8 

5.4 

(1863) 

8 

4.9 

(1676) 

8 

4.4 

(1520) 

8 

9.6 

9.5 

9.4 

9.1 

8.7 

8.2 

7.6 

(3322) 

(3244) 

(3125) 

(2903) 

(2696) 

(2488) 

(2285) 

L 

L 

L 

L 

L 

L 

L 

11.9 

12.1 

12.3 

12.3 

12.2 

12.0 

11.8 

(3949) 

(3880) 

(3756) 

(3497) 

(3289) 

(3118) 

(2965) 

U”-*-  MM 

l 

L 

l 

L 

L 

L 

L 

1. 

15  4 

16.4 

17.7 

19.4 

20.2 

20.4 

20.4 

1.  1 

(4347) 

(4324) 

(4246) 

(4045) 

(3829) 

(3625) 

(3451) 

l 

i 

L 

L 

L 

L 

l 

■.v 

20  8 

22.8 

25.8 

30.3 

32.5 

33.6 

34.0 

O 

□0 

(4428) 

(4467) 

(4339) 

(4120) 

(3902) 

(3703) 

'  ^ 

L 

L 

L 

L 

L 

L 

L 

31.3 

35  4 

41.8 

52.9 

57.8 

60.2 

61  2 

(4550) 

(4568) 

(4610) 

(4649) 

(4389) 

(4137) 

(3908) 

l 

L 

L 

L 

l 

L 

L 

63.0 

73  9 

91.2 

122.9 

137  9 

143.0 

145.2 

(4870) 

(4884) 

(4922) 

(5040) 

(4775) 

(443G) 

(4150) 

l. 

l 

L 

L 

l 

L 

L 

c.  Fuel:  HMX. 


liquid  oxide,  8  *  T1O2  (s). 


7.0 

6.8 

6.4 

5,7 

49 

4  3 

3.6 

(2492) 

(2396) 

(2235) 

(1977) 

(1706) 

(1481) 

(1290) 

•V  \‘:-i 

L 

L 

L 

e 

8 

S 

8 

.  - 

96 

9.5 

9.3 

8.8 

8  1 

7  3 

6.5 

r  *■ 

(3311) 

(3222) 

(3084) 

(2808) 

(2519) 

(2218) 

(19561 

w*.  ■ 

L 

l 

l 

L 

L 

L 

8 

11.9 

12.1 

12.2 

12  1 

11  8 

11  4 

10.7 

(3944) 

(3049) 

(3731) 

(3431) 

(3175) 

(2943) 

(2697) 

•  ••  V 

L 

L 

L 

L 

l 

L 

L 

, 

15.4 

16  3 

17  6 

19.1 

19.6 

19.5 

19.2 

(4340) 

(4313) 

(4225) 

(3999) 

(3744) 

(3491) 

(3267) 

1 

L 

L 

L 

L 

L 

l 

L 

20  8 

22.7 

25  7 

29.9 

31  7 

32  2 

32  0 

*"»/// .V 

(4396) 

(4406) 

(4429) 

(4287) 

(4032) 

(3765) 

(3514) 

L 

L 

l 

l 

l 

l 

L 

31.2 

35.2 

41.4 

52.1 

56.2 

57  5 

57.4 

‘  *  *  *  *  \ 

(4536) 

(4540) 

(4560) 

(4579) 

(4285) 

(3979) 

(3691) 

* — , 

L 

L 

l 

L 

L 

L 

l 

62  8 

73.3 

89.9 

119.9 

133.7 

135  7 

1350 

,  .  ■  .■ 

(4852) 

(4849) 

(4857) 

(4910) 

(4639) 

(4236! 

(3879) 

-''i 

L 

L 

l 

L 

L 

L 

L 

11 


.  1.*  *. 


NWCTP6544 


TABLE  5.  (Contd.). 


C  +  M 

C/M 

0.1 

02 

0.4 

1.0 

2.0 

4.0 

10.0 

d.  Fuel:  Perrtolite. 


.1 

7.1 

6.9 

66 

6.1 

5.5 

5.0 

(2510) 

(2431) 

(2301) 

(2102) 

(1881) 

(1698) 

L 

l 

L 

B 

8 

B 

■2 

9.6 

9.5 

9  4 

9.1 

8.7 

8.3 

(3324) 

(3247) 

(3130) 

(2914) 

(2715) 

(2517) 

L 

L 

L 

L 

L 

L 

.4 

11.9 

12.1 

12.3 

12.3 

12.2 

12.1 

(3951) 

(3883) 

(3762) 

(3505) 

(3297) 

(3130) 

L 

L 

L 

l 

L 

L 

.0 

15.4 

16.4 

i  7.6 

19.4 

20.1 

20.3 

(4334) 

(4313) 

(4236) 

(4037) 

(3822) 

(3613) 

L 

L 

l 

L 

L 

L 

:.o 

20.7 

22.6 

25.6 

30.3 

32.5 

33.4 

(4385) 

(4384) 

(4389) 

(4311) 

(4091) 

(3863) 

l 

L 

L 

L 

l 

L 

1.0 

31.1 

35.0 

41.1 

57.8 

59.8 

(4523) 

(4514) 

(4508) 

(4517) 

(4343) 

(4071) 

L 

L 

L 

L 

L 

L 

0.0 

62.3 

72.4 

88.3 

117.3 

137.8 

142.1 

(4839) 

(4822) 

(4800) 

(4778) 

(4710) 

(4345) 

L 

L 

L 

L 

L 

L 

e.  Fuel: 

Comp  B. 

0.1 

7.1 

69 

65 

6.0 

5.3 

4.7 

(2503) 

(2418) 

(2277) 

(2056) 

(1817) 

(1618) 

L 

L 

l 

B 

B 

B 

0.2 

9.6 

9.5 

9.4 

9.0 

8.5 

7.9 

(3319) 

(3238) 

(3113) 

(2876) 

(2648) 

(2414) 

l 

L 

L 

L 

L 

L 

04 

119 

12.1 

12.3 

12.3 

12.1 

11.8 

(3948) 

(3877) 

(3750) 

(3477) 

(3252) 

(3065) 

l 

l 

L 

L 

l 

l 

1.0 

15.5 

16.4 

17.7 

19  3 

20.0 

20.0 

(4326) 

(4302) 

(4220) 

(4009) 

(3776) 

(3546) 

L 

L 

L 

l 

L 

l 

2.0 

20  8 

22.7 

25.6 

30  3 

32.2 

32.8 

(4373) 

(4362) 

(4350) 

(4270) 

(4023) 

(3771) 

l 

l 

L 

l 

l 

l 

4.0 

31.1 

35.0 

40.9 

51.4 

57.4 

58.6 

(4S09) 

(4488) 

(4460) 

(4421) 

(4264) 

(3949) 

L 

l 

L 

L 

L 

L 

100 

62  2 

72,2 

87.5 

114.9 

135.4 

139.0 

(4822) 

(4791) 

(4743) 

(4655) 

(4567) 

(4187) 

l 

L 

L 

l 

l 

L 

0545) 

B 

7.7 

(2323) 

L 

11.8 

(2983) 

L 

20.3 

(3432) 

L 

336 

(3647) 

L 

603 

(3811) 

L 

142.1 

(4000) 

L 


4.2 

(1452) 

B 

73 

(2188) 

L 

11.5 

(2893) 

L 


19.9 

(3344) 

l 

32.7 

(3520) 

L 

58.2 

(3632) 

L 

135.6 

(3756) 

L 


L  *  liquid  oxide,  B  »  riO^O) 
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TABLE  5.  (ContcL). 


cm 

0.1 

0.2 

J  04 

1.0 

20 

4.0 

10.0 

f.  Fuel:  TNT. 


0.1 

7.1 

7.0 

6.7 

6.2 

5.9 

5.5 

5.1 

<2524) 

(2459) 

(2352) 

(2147) 

(2013) 

0862) 

(1736) 

L 

L 

L 

L 

B 

B 

3 

0.2 

9.6 

9.6 

9.5 

9.3 

9.1 

8.9 

8.6 

(3334) 

(3265) 

(3163) 

(2982) 

(2832) 

(2695) 

(2565) 

L 

L 

L 

L 

L 

L 

L 

0.4 

12.0 

12.2 

12.4 

12.6 

12.5 

12.4 

12.3 

(3954) 

(3891) 

(3780) 

(3547) 

(3355) 

(3205) 

(3082) 

l 

L 

L 

L 

L 

L 

L 

1.0 

15  5 

16.4 

17.8 

19.6 

20.4 

20  5 

20.2 

(429S) 

(4264) 

(4195) 

(3995) 

(3776) 

(3537) 

(3292) 

L 

L 

L 

L 

L 

L 

L 

2.0 

20.7 

22  5 

25.2 

29  5 

32.3 

33.3 

31.8 

(4331) 

(4281) 

(4201) 

(4037) 

(3874) 

(3641) 

(3252) 

L 

L 

L 

L 

L 

L 

L 

4.0 

30  8 

34  3 

39  4 

48.1 

53  6 

56.1 

54  8 

(4463) 

(4404) 

(4352) 

(4194) 

(3979) 

(3682) 

(3259) 

L 

CL 

CL 

CL 

CL 

CL 

CL 

10.0 

60.7 

69.2 

82  7 

106  0 

120  7 

126.2 

122.3 

(4812) 

(4791) 

(4732) 

(4531) 

(4244) 

(3633) 

(3286) 

CL 

CL 

CL 

Cl 

CL 

CL 

CL 

g.  Fuel:  Hydrazine. 


0.1 

7  2 

7  1 

6.9 

66 

C  ’ 

5.8 

54 

(2510) 

(2437) 

(2320) 

(2143) 

(1961) 

(1809) 

(1684) 

L 

L 

L 

B 

B 

6 

B 

0.2 

9  7 

9.8 

9  9 

10.0 

99 

9.7 

94 

(3303) 

(3227) 

(3118) 

(2921) 

(2  730) 

(2589) 

(2436) 

l 

L 

L 

l 

L 

1. 

L 

0  4 

12  3 

12  7 

13  3 

140 

14  2 

14  0 

13  5 

(3899) 

(3805) 

(3659) 

(3386) 

13141) 

(2895) 

(2644) 

L 

L 

L 

L 

L 

L 

L 

1.0 

17.0 

18.8 

21.2 

23.4 

22.8 

20.8 

18  3 

(4277) 

(4146) 

(3983) 

(3555) 

(3005) 

(2457) 

(1979) 

. 

L 

L 

L 

l 

L 

L 

A 

2.0 

23  4 

27  0 

32  0 

39  0 

39  1 

29  9 

28.2 

(4252) 

(4152) 

(3991) 

(3626) 

13042) 

(2047) 

(1739) 

L 

L 

L 

L 

L 

A 

AB 

4.0 

36  2 

43.3 

53.4 

68.0 

70.9 

55  3 

43.7 

(4373) 

(4248) 

(4104) 

(3735) 

(3101) 

(2062) 

(1443) 

L 

L 

NL 

NL 

NL 

NL 

AO 

10.0 

75  6 

93  8 

119  9 

157  9 

165  0 

120.3 

93.1 

(4650) 

(4482) 

(4374) 

(3927) 

(3148) 

(2068) 

(1291) 

L 

L 

NL 

NL 

NL 

NL 

AN 

L  «  liquid  oxide,  A  =  T13C5  (s),  B  T1O2  (s),  C  -  TiC,  N  =  TiN 
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TABLE  5.  (Contd.). 


C  +M 

C/Af 

0.1 

0.2 

0,4 

1.0 

2.0 

40 

10.0 

h. 

Fuel:  Eth) 

(lene  oxide 

0.1 

7.3 

7.4 

7.5 

7.8 

7.9 

7.9 

(2580) 

(2568) 

(2551) 

(2502) 

(2487) 

(2475) 

t 

L 

L 

L 

L 

L 

0.2 

9  8 

10.0 

10.2 

10.9 

11  1 

11.2 

(3304) 

(3314) 

(3261) 

(3103) 

(3034) 

(2963) 

L 

L 

L 

L 

L 

L 

L 

0.4 

12.3 

12.8 

13.4 

14.3 

14.5 

14.2 

13.7 

(3929) 

(3857) 

(3743) 

(3502) 

(3246) 

(3975) 

(2711) 

l 

L 

L 

L 

L 

L 

L 

1,0 

16  6 

10.2 

20.3 

22.7 

23.5 

22.4 

18  9 

(4194) 

(4083) 

(3903) 

(3551) 

(3220) 

(2760) 

(2164) 

l 

l 

L 

CL 

CLN 

CN 

CG 

2.0 

22.6 

26.4 

29.3 

35.0 

37  3 

32? 

29  6 

(4207) 

(4085) 

(3984) 

(3676) 

(3252) 

(2578) 

(2133) 

l 

CL 

CL 

CL 

CL 

CG 

CG 

4.0 

34.1 

39.6 

47.9 

59.9 

62  6 

57  0 

51.1 

(4374) 

(4318) 

(4204) 

(3820) 

(3268) 

(2612) 

(2115) 

CL 

CL 

CL 

CL 

C 

CG 

CG 

100 

68.9 

83.3 

105  1 

136.2 

1368 

1298 

115  8 

(4741) 

(4676) 

(4535) 

(4022) 

(3259) 

(2648) 

(2107) 

1 _ Et _ 

_£L_J 

CL 

CL 

_ £ _ 

CG 

CG 

i.  Fuel:  Carbon. 

0 1 

7.4 

75 

7  6 

78 

80 

8.1 

(2616) 

(2627) 

(2643) 

(2567) 

(2678) 

(2671) 

L 

L 

L 

L 

L 

L 

02 

9.8 

99 

10.0 

10  2 

90 

78 

(3404) 

(3386) 

(3350) 

(3241) 

(2657) 

(2260) 

L 

L 

L 

L 

L 

L 

L 

0.4 

12.1 

12  3 

12  6 

12  1 

10  5 

8  4 

70 

(3979) 

(3929) 

(3792) 

(3336) 

(2821) 

(2305) 

(1960) 

L 

L 

L 

N 

CG 

CG 

GN 

1 0 

15.2 

15.3 

15.2 

13  8 

11  2 

80 

58 

(4067) 

(3908) 

(3804) 

(3588) 

(3010) 

(2221) 

(1709) 

L 

CL 

CL 

CG 

CG 

CG 

GN 

2.0 

19.2 

19  1 

20.0 

150 

11.7 

78 

5  1 

(4168) 

(4126) 

(4473) 

(3788) 

(3128) 

(2170) 

(1558) 

CL 

CL 

C 

CG 

CG 

CG 

GN 

4.0 

26  5 

26  1 

24.5 

16.4 

12  1 

7  7 

3.5 

(4414) 

(4375) 

(4716) 

(3962) 

(3212) 

(2135) 

(1421) 

C.L 

CL 

C 

CG 

CG 

CG 

AGN 

10.0 

48.2 

47.8 

36.4 

189 

12  3 

76 

2.4 

(4792) 

(5133) 

(5014) 

(4165) 

(32  76) 

(2110) 

(1365) 

CL 

C 

C 

CG 

CG 

C.G 

AGN 

L  «  liquid  oxide.  C  »  TiC.  G  «  graplute.  N  »  Tin,  A  ■  11305(1) 
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TABLE  5.  (Contd.). 


cm 

C  r  M 

0.1 

0-2  . 

. 

0.4 

1.0 

4.0 

10.0 

j.  Fuel: 

Hexane. 

0.1 

7,6 

7.8 

8.2 

8.9 

9.3 

9.6 

9.6 

(2650) 

(2691) 

(2749) 

(2824) 

(2850) 

(2833) 

(2771) 

L 

L 

L 

L 

L 

L 

L 

0.2 

10  1 

104 

10.9 

11.4 

11.0 

10  1 

93 

(3387) 

(3373) 

(3336) 

(3127) 

(27571 

(2391) 

(2082) 

t 

L 

L 

L 

L 

L 

L 

04 

12.7 

134 

144 

147 

12.7 

108 

9.1 

(3902) 

(3794) 

(3572) 

(3037) 

(2359) 

(1925) 

OS52) 

L 

L 

L 

IN 

CN 

GN 

GN 

'.0 

17  4 

18  9 

20.8 

20.2 

17  7 

15  3 

12  6 

(4005) 

(3805) 

(3629) 

(2826) 

(2184) 

(1760) 

(1377) 

L 

CIN 

CLN 

C 

CG 

GN 

A8G 

20 

23  6 

26.6 

307 

30.5 

27  3 

23  3 

19.2 

(4075) 

(3972) 

(374Q) 

(2891) 

(2162) 

(1674) 

(1297) 

CL 

CL 

CL 

CG 

CG 

GN 

BG 

40 

35.8 

42  3 

52.2 

51.3 

46  5 

39  4 

31.9 

(4293) 

(4177) 

(4042) 

(2951) 

(2151) 

(1622) 

(1205) 

Cl 

CL 

C 

CG 

CG 

GN 

BG 

10.0 

73  4 

90  2 

111  5 

113.7 

104.3 

88  2 

64.2 

(4623) 

(4481) 

(4276) 

(3009) 

(2146) 

(1586) 

(1025) 

CL 

CL 

C 

CG 

CG 

GN 

UGN 

L  o  liquid  oxidu.  A  «  T13Q5  (s),  B  ••  Ti-jO*j(s),  C  »  TiC,  G  »  graphitu.  N  »  TiN 


In  Tublo  6  is  given  the  sul  of  product  yields  for  TNT-titunium-air  ut  CIM  =  1 .0,  chosen 
us  representative.  Mole-percent  in  the  vupor  is  given  for  eudi  species  that  reaches  at  leust 
0.1%  at  some  point.  Numbers  of  moles  of  condensed  phases  are  also  given.  As  with  Table  4, 
the  composition  of  the  liquid  oxide  phase  is  given  aH  x  in  the  empirical  formula  TiOx. 

Product-yield  data  are  available  for  all  the  systems  run  but  have  not  been  included  in 
this  report  due  to  the  bulk  ofduta  involved.  They  are  uvuiluble  from  the  author. 


DISCUSSION 

Figure  1  shows  the  effect  of  total  concentration  (C  4  M)  on  udiabatic  lemporuturo  for 
pure  titanium  in  air  as  well  as  for  two  represontutive  C/M  ratios  for  the  oxygon-deficient 
hexane  and  the  oxygen-rich  nitrocellulose. 
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TABLES.  Combustion  of  TNT-Titanium  in  Air  at  C1M  -  1.0: 

Overpressure,  Adiabatic  Temperature,  and  Product  Yields. 


Property  / 
system 

C  +  M,  kg/m3 

0.1 

0.2 

0.4 

1.0 

2  0 

4.0 

10.0 

Overpressure,  bars 

5.2 

9.3 

12.5 

19.6 

29.5 

48.1 

106.0 

Temperatve,  K 

2147 

2982 

3547 

3995 

4037 

4194 

4531 

Mole  %: 

Ar 

0.99 

0.97 

0.88 

0.65 

0.44 

0.29 

0.14 

Ti 

0.03 

4.47 

6.80 

7.84 

TiO 

0,31 

4.35 

3.22 

2.41 

2.57 

Ti02 

0.02 

0.02 

TI  + 

001 

0.01 

CO 

0.01 

1,82 

,979 

2420 

33  46 

37.91 

40.95 

co2 

3.78 

5  58 

3.62 

0.62 

H 

0.14 

1.53 

7.02 

RS5 

10.01 

12.30 

OH 

008 

1.02 

2.19 

0.91 

001 

0.01 

Hj 

0.09 

074 

4.03 

7.65 

10.06 

12.22 

H20 

1  31 

1.97 

'  2.19 

0.87 

001 

001 

001 

NO 

0.99 

308 

3.14 

0.62 

n2 

77.86 

75  58 

69  80 

5573 

41  72 

31.12 

21.69 

O 

003 

098 

2  66 

0  89 

°2 

14.94 

878 

.3.13 

0.08 

CN- 

0.01 

0.01 

CN 

0  09 

0.25 

0  46 

c2h 

001 

0.06 

0.20 

CjN 

0.02 

006 

HCN 

0  34 

0.97 

1  52 

HNCO 

001 

HCO 

001 

001 

0.03 

0  06 

QH2 

003 

0  11 

Liquid  oxide. 

Moles 

1  04 

1.70 

2.58 

6  48 

13  88 

22  60 

45.39 

x  in  TiO* 

2.00 

1  94 

1.69 

1  18 

072 

0  57 

0  56 

Moles  TiC  (1) 

6  17 

29.86 

The  behavior  of  the  pure  metal  is  quite  unlike  that  previously  seen  for  aluminum 
(where  a  maximum  in  T  versus  C  +  M  occurred  at  approximately  the  AI2O3  stoichiometry; 
see  Ref.  4),  With  titanium  there  is  a  steady  increase  in  temperature  with  concentration,  the 
result  of  the  existence  of  a  number  of  oxides,  the  formation  of  each  being  an  exothermic 
process.  Again,  wherous  aluminum  showed  formation  of  AIN  at  high  concentrations,  TiN 
does  not  form  here  as  a  product  owing  to  the  lesser  high-temperature  stability  of  TiN.  The 
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curves  for  the  explosive  (relatively  oxygen-rich)  fuels  are  very  similar  to  those  for  the  metal 
alone. 

With  the  oxygen-deficient  fuels,  first  TiN  and  then  TiC  and  graphite  appear  as 
products  in  the  carbon-rich  region.  These  tend  to  cause  a  falling  off  of  temperature,  as  is  seen 
in  the  case  of  hexane  in  Figure  1 . 

Figure  2  shows  the  strong  effect  of  total  concentration  on  overpressure.  Nitrocellulose 
has  been  chosen  as  representative;  all  the  fuels  (except  carbon)  show  a  similar  pattern.  The 
large  increase  in  number  of  moles  of  gaseous  products  combined  with  a  much  less  profound 
change  in  temperature  (Figure  1)  is  responsible  for  the  monotonic  increase  in  pressure. 

The  effect  of  C/M  on  temperature  is  shown  in  Figure  3.  In  general,  titanium  has  a 
higher  heating  value  than  any  of  the  fuels  used;  hence,  there  is  a  generally  observed  decrease 
in  temperature  aB  titanium  is  replaced  by  fuel.  At  high  concentrations  with  oxygen-rich 
fuels,  the  oxygen  supplied  by  the  fuel  causes  the  temperature  to  rise  at  first,  resulting  in  a 
shallow  maximum  in  the  curve.  The  considerable  break  downward  in  the  curves  for  hexane 
again  corresponds  to  the  production  of  TiC  and  graphite. 

Commonly,  overpressure  increases  with  C/M,  as  seen  in  Figure  4,  owing  to  the 
increase  in  number  of  moles  of  gases  (CO  +  H2)  as  titanium  is  replaced  by  fuel.  The 
maximum  seen  with  hexane  reflects  the  sharp  drop  in  temperature,  referred  to  above. 

The  nature  of  the  fuel  has  little  effect  on  the  adiabatic  temperature  for  the  oxygen- 
rich  fuels.  For  the  oxygen-deficient  fuels  there  is  a  general  increase  in  temperature  as  the 
oxygen  balance  increases  (becomes  less  negative),  as  is  shown  in  Figure  5,  presumably  due  to 
relieving  the  oxygen  lack.  Pressure  is  surprisingly  insensitive  to  the  nature  of  the  fuel, 
except  for  the  case  of  carbon  which  shows  higher  temperatures  and  markedly  lowet  pressures 
as  a  result  of  a  much  smaller  quantity  of  gaseous  product  thun  with  the  hydrogen-containing 
fuels.  Carbon  also  shows  a  very  different  set  of  products  from  those  represented  in  Table  6. 
With  carbon  at  high  concentrations  are  seen  large  amounts  of  CN,  C2N  and  C3  and  even 
several  tenths  mole-percent  of  the  ions  as  a  result  of  the  high  temperatures  experienced, 
combined  with  the  absence  of  hydrogen. 


APPROXIMATIONS 

Athow  has  examined*1  the  ideal-gas  approximation  for  computation  on  internal 
explosions  in  the  presence  of  the  reactive  metal  aluminum.  In  those  systems  it  appears  that 
the  only  gas  below  its  critical  temperature  is  aluminum  vapor;  and  for  aluminum7  the 
reduced  pressure  can  be  seen  not  to  exceed  10-4  and  the  reduced  volume  not  to  be  less  than 
200  so  that  no  appreciable  deviation  from  ideality  should  occur.  It  is  anticipated  that  the 
same  conclusions  can  be  made  regarding  the  titunium  systems,  although  no  estimates  of  the 
critical  properties  have  been  made  and  thus  no  estimates  of  errors  attempted. 


1*  L.  K.  Athow,  "Real  Gas  Considerations  for  Determining  Physical  und  Thurmodynumic  Properties  of 
Guses  Involved  in  the  Prediction  of  the  Effects  of  internal  Explosions."  M.  S.  Thesis,  Naval  Postgraduate  School,, 
Monterey, Calif.,  June  1982. 
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It  was  pointed  out  by  Smith18  that  radiation  effects  should  be  negligible  when 
detonation  takes  place.  When  the  combustion  is  a  deflagration,  however,  the  adiabatic 
assumption  is  expected  to  be  less  reliable.  A  comparison  of  the  experimental  data  on  dust 
explosions  reported  by  the  Bureau  of  MineslS  would  seem  to  show  that  the  computed  results 
from  the  present  report  may  be  high  by  a  factor  of  about  two.  It  has,  however,  been  pointed 
out  by  Baker  et  alM  that  the  Hartmann  apparatus  used  by  the  Bureau  of  Mines  investigators 
seriously  underestimates  the  maximum  overpressures  to  be  expected.  On  this  basis  it  is 
anticipated  that  the  adiabatic  results  will  show  much  less  extreme  negative  departures  from 
realistic  maximum  overpressures. 


13  D.  E.  Smith,  "Attenuation  Effect*  of  Thermal  Radiation  on  Internal  Blast  Overpressure."  M.  S.  Thesis. 
Naval  Postgraduate  School,  Monterey ,  Calif.,  December  1979. 

13  Bureau  of  Minea.  Explatibility  of  Mtial  PounUrt,  by  M.  Jacobson,  A.  R.  Cooper  and  J.  Nagy.  U.S. 
Department  of  Interior,  Washington,  D.C.,  1964.  (R1 6516,  publication  UNCLASSIFIED.) 

1*  W.  E.  Baker,  P.  A.  Cox,  P.  S.  Westine,  J.  J.  Kulesz,  and  R.  A.  Strehlow.  Expiation  Hazard)  and 
Evaluation.  Amsterdam,  Elsevier,  1982,  pp.  260-6 1 . 
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Appendix  A 

DESCRIPTION  OF  THE  COMPUTER  PROGRAM 


The  program  itself,  which  is  listed  in  Appendix  B,  is  written  for  the  HP  9845A 
computer.  Various  sections  are  referred  to  by  their  labels. 


MAIN  PROGRAM 
Input  Section 

Fuel  (accessible  from  line  1900).  Enter  formula  and  internal  energy  of  formation. 
Computes  formula  mass;  allows  for  zero  C  or  H  or  no  fuel.  Resets  flags. 

Cone  (accessible  from  line  1870).  Resets  counters.  Enter  C  +  M  (total  concentration 
in  kg/m3)  and  then  C/M,  as  asked  for. 

Temp  (for  entering  total  temperature  manually).  May  be  accessed  by  use  of  special 
function  key  k4.  After  at  least  two  trials,  Temp  may  be  bypassed  and  interpolation  used  to 
find  the  new  temperature.  Pressing  k5,  so  that  At$  a  "Y”,  after  two  iterations  allows  for 
automatic  interpolation.  If  needed,  the  automatic  interpolation  may  be  stopped  by  pressing 
PAUSE,  then  k5  (which  returns  At$  to  "N”),  and  finally  k4  for  manual  setting  of 
temperature. 


Computation  Section 

Calc  calls  up  the  computational  subroutines  Eq,  Tical,  and  Ex. 

Eq  is  a  subroutine  of  the  main  program  to  evaluate  equilibrium  constants  of 
formation  of  each  of  the  40  chemical  species.  Kp  is  first  computed  from  the  stored  parameters 
and  then  converted  to  Kn  (expressed  in  mole  numbers).  Kp  or  Kn  is  defined  as  the  ratio  of  the 
activity  of  the  species  to  the  product  of  the  activities  of  C,  Ti  vapor,  H2,  N2,  and  O2,  each  raised 
to  the  appropriate  power.  For  and  T13O5,  >/Kn  is  computed  so  as  to  avoid  overflow  at 
lower  temperatures.  For  condensed  phases,  the  standard  state  is  the  pure  phase;  for  Kn,  the 
standard  state  for  gases  is  1  mole  (in  1  m3). 


Ileal  is  the  master  subroutine  which  carries  out  equilibrium  and  energy  calculations. 
Results  are  displayed  as  "dU”  (net)  and  "T  high”  or  "T  low”;  dU  =  0  is  desired  for  convergence. 
A  new  temperature  approximation  is  performed  automatically  by  interpolation,  based  on  the 
previous  T  and  dU  values;  or  else  it  is  entered  manually  with  k4.  A  more  detailed  description 
of  Tical  will  be  found  below. 
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Output 

At  Ez,  results  of  the  calculations  are  printed:  temperature,  overpressure,  mole 
numbers  of  all  products  and  the  diagnostic  features  described  under  DIAGNOSTICS,  below. 
Options  are  allowed  at  this  point  for  changing  C/M,  C  +  M  or  for  a  new  fuel. 


SUBROUTINE  TICAL 

Tical  is  the  major  computational  subroutine  whose  task  la  to  find  the  numbers  of 
moles  of  the  products  present  at  equilibrium  at  the  selected  temperature.  The  conditions  to  be 
satisfied  (other  than  for  the  trivial  case  of  argon)  are  the  atom  balances  in  C,  H,  N,  O,  and  Ti 
and  the  establishment  of  chemical  equilibrium  between  each  compound  and  its  component 
elements  in  their  reference  states. 

The  master  variables  (all  in  mole  numbers)  are  X  =  Voj,  Y  =  VH2.  Z  -  VN^ 
Acc  =  activity  of  carbon  (standard  state  -  graphite),  and  Ti  =  Ti  metal  vapor.  Of  those,  Y  is 
always  computed  in  closed  form;  from  one  to  four  of  the  remaining  are  found  as  unknown 
parameters  in  the  subroutine  Newt  by  using  the  Newton-Kaphson  method.  The  actual 
number  of  unknowns  is  equal  to  four,  reduced  by  the  number  of  condensed  phases  present. 
Possible  condensed  phases  are:  Ti  (solid),  TiO  (solid),  Ti203  (solid),  T^Ob  (solid),  TiOg  (solid), 
liquid  oxide  solution  (designated  Lox  in  the  program),  C  (solid),  TiC  (solid  or  liquid),  TiN 
(solid  or  liquid).  The  presence  of  each  condensed  phase  is  indicated  by  a  flag,  to  be  set  as 
described  under  Flags. 

Based  on  the  values  of  the  master  variables  and  the  equilibrium  constant  of 
1  formation,  the  mole  number  of  each  species  is  computed.  Then  the  material  balance  in  the 

elements  O,  N,  0,  and  Ti  is  written  in  terms  of  these  mole  numbers.  When  liquid  oxide  is 
present,  the  stoichiometric  condition  is  that  the  sum  of  the  activities  (that  is,  mole-fractions) 
of  the  components  of  the  solution  should  add  to  unity;  this  condition  replaces  the  atom  balance 
in  titanium  for  this  case.  There  thus  results  a  set  of  up  to  four  simultaneous  nonlinear 
equations;  this  set  is  the  basis  of  the  Newton-Raphson  scheme  to  find  the  unknown 
L  parameters. 

At  the  conclusion  of  an  iteration,  the  newly  generated  values  of  the  master  variables 
are  used  to  repeat  the  calculations.  In  favorable  situations,  each  iteration  results  in 
improvement  (although  temporary  divergence  sometimes  occurs).  Iteration  is  repeated  until 
the  stoichiometric  errors  are  less  than  one  part  in  ten  thousand. 

Newt  solves  the  set  of  simultaneous  nonlinear  equations  needed  to  find  the  master 
variables.  The  Newton-Raphson  method  is  used.  In  this  subroutine  it  is  necessary  to  find  a 
number  of  derivatives  numerically  by  observing  the  effect  of  a  fractional  change  in  each 
variable.  This  fractional  change  is  set  initially  at  0.1.  It  is  found  that  divergences  which 
would  otherwise  occur  when  the  errors  are  large  can  be  averted  by  altering  this  fractional 
change  to  0.5.  This  is  accomplished  by  the  use  of  special  function  key  k8,  when  the  computer 
is  in  a  pause  mode  or  is  waiting  for  input.  Depressing  k8  a  second  time  will  change  the 
fraction  back  to  0. 1 .  The  fractional  errors  in  the  stoichiometric  conditions  may  be  observed  by 
the  use  of  kl  (TRACE  VARIABLES  Yn(*)/EXECUTE)  after  which  will  be  displayed  Y(l), 
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Y(2),  etc.,  which  are  these  fractional  errors,  The  display  of  these  will  remain  on  the  screen  if 
PUT  ALL  is  locked  down. 

Approx:  To  begin  the  calculation  an  initial  approximation  of  the  master  variables  is 
needed.  In  this  approximation  an  arbitrary  hierarchy  of  oxygen  and  nitrogen  uptake  is 
assumed.  Oxygen  is  assumed  to  be  taken  up  in  the  order  CO,  titanium  oxides,  H2O,  CO2.  If 
there  is  insufficient  oxygen  to  form  TiO,  then  TiN  is  assumed  to  be  present.  With  excess  C, 
C  (solid)  or  TiC  are  assumed. 

It  has  been  found  that  Approx  gives  quite  satisfactory  values  of  the  master  variables 
at  temperatures  below  about  3000  K;  at  higher  temperatures  the  initial  approximation  may 
be  in  considerable  error.  After  each  entry  into  Temp  the  question  is  asked,  "Do  you  want 
'Approx'?".  The  default  condition  (obtained  by  pressing  CONT)  is  "Y"  for  a  new  concentration 
(C/M  or  C  +  M)  and  "N"  otherwise. 

Flags  sets  the  flags  for  condensed  phases  and  gives  the  values  of  the  temporary 
variables  used  in  Newt.  The  criteria  for  the  presence  of  a  condensed  phase  are  that  the 
quantity  of  the  phase,  if  previously  computed,  be  nonnegative  and  that  the  formation  constant 
be  satisfied  or  exceeded.  In  the  case  of  liquid  oxide,  the  formation-constant  requirement  is 
replaced  by  the  requirement  that  the  temperature  be  above  the  melting-point  curve  (T  = 
2143  is  used  for  simplicity;  see  text  for  discussion  of  this  point)  and  that  the  sum  ol  the 
activities  of  the  components  of  the  solution  be  no  less  than  unity. 

Two  important  indexes  are  set  by  Flags.  liflag  gives  the  number  of  metal-plus-oxide 
phases  (no  greater  than  two);  for  liquid  oxide  solution,  liflag  =  1  except  that,  when  TiC, 
C  (solid),  and  liquid  oxide  are  all  present,  liflag  =  0.  The  index  li  gives  the  number  of 
unknowns  to  be  sought  in  Newt.  Finally  Flags  gives  initial  values  of  the  variables  to  be  used 
in  Newt. 

Fx  gives  the  fitting  functions  for  Newt.  There  are  three  branches,  depending  on  the 
value  of  liflag.  Certain  of  the  master  variables  are  computed  in  Fx.  After  Dlff  is  called,  the 
current  mole  numbers  of  all  species  are  generated,  The  subroutine  returns  to  Newt  the 
variuble  Fx,  which  is  the  fractional  error  in  the  stoichiometry  in  whichever  element  Newt  is 
considering  at  the  time. 

Diff  first  computes  those  master  variables  which  were  not  found  in  Fx.  Then  Spec 
and  Oxides  are  called  to  compute  the  mole  numbers  of  all  species,  and  finally  the  errors  in 
stoichiometry  for  all  elements  (and  the  sura  of  activities  when  liquid  oxide  is  present)  are 
computed. 

Spec  computes  the  mole  numbers  of  all  gaseous  species,  given  the  current  values  of  X, 
Z,  Acc,  and  Ti.  Y  =  VH2  is  computed  in  closed  form  in  Spec;  it  is  needed  for  computations  on 
hydrogen-containing  species.  For  each  species  the  appropriate  formation  equilibrium 
constant  is  used. 

Oxides  computes  the  mole  numbers  of  condensed  metal  and  oxide  phases.  For  two 
such  phases  (liflag  =  2)  the  atom  balance  conditions  for  titanium  and  oxygen  are  used;  for 
liflag  =  1,  only  the  titanium  balance  is  used;  for  liflag  =  0  the  oxygen  balance  is  used. 
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Nwi  serves  to  route  the  computation  to  the  proper  portions  of  Tical.  On  first  entry, 
the  flags  will  have  been  set  and  control  is  sent  to  Newt;  except,  if  Ii  =  0,  provision  exists  for  a 
closed-form  solution.  In  either  case,  Flags  is  called  again.  If  any  flag  has  changed,  the 
computation  must  be  repeated  for  the  new  set  of  condensed  phases;  othorwise  preparation  is 
made  for  exit  from  Tical  at  the  line  Energy , 

Energy  computes  from  the  stared  parameters  the  molar  internal  energy  of  each 
component  and  then  tho  total  change  in  internal  energy  (dll)  from  tho  starting  materials, 
taking  into  account  tho  mole  numbers  of  each  species  present.  If  the  system  is  at  the  molting 
point  of  any  of  the  condensed  phu  jos  present,  tho  amounts  of  solid  und  liquid  are  computed 
from  the  energy  balance,  using  the  known  energy  of  fusion.  Control  is  then  returned  to  the 
main  program. 


diagnostics 

Suml  is  culled  at  the  end  of  ouch  run  to  report  the  following  itoms;  (1)  a  chock  on  the 
material  balance  of  ouch  element;  (2)  the  activity  of  carbon  (Acc)  und  tho  sum  of  tho  activities 
of  the  components  of  the  liquid  oxide  solution;  (3)  a  comparison  of  the  computed  amounts  with 
tho  equilibrium  constants  of  formation  for  each  condonsed  phase;  and  (4)  the  activities  (mole* 
fractions)  of  ouch  component  in  tho  liquid  oxide  solution. 

Sum  is  called  whenever  special  function  key  kO  has  been  depressed  onco.  At  each 
emergence  from  Newt  tho  mole  numbers  of  ail  species  are  presented  in  condensed  form;  this  is 
followed  by  Suml,  the  output  of  which  has  been  just  described.  These  chucks  are  of  particular 
interest  in  troubleshooting.  Depressing  kO  a  second  time  will  cancel  calling  up  Sum. 

Depressing  kl  will  cuuse  execution  of  TRACE  VARIABLES  Yn(*).  Then  during  each 
iteration  in  Newt  the  rotative  error  functions  used  to  tost  convergence  will  be  displayed. 
Euch  Yn(*)  must  drop  below  0,0001  for  convergence.  This  feuturu  muy  bo  turned  off  by 
executing  NORMAL. 
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Appendix  B 

PROGRAM  LISTING  FOR  "TIF4H 


It 

9 

V, 

tt 


FILE  MANE  'TIF4'  22  Fob  84  i  Tl  ♦  air  ♦  fill)  cintlniiit  Uqild  phatu 
Adiabatic  iiirpritiari  calcilatlin  (cinaiant 


l\ 


IWIjtJ  TO  *TMT1 


a. 


0 

111 

N 

II 

I 

I 

171 

1BI 

170 

211 

g| 

231 


to 


DEM  M.l 

INTEKR  C(SiV. X.i.K , li.N«,CijNMc  t . In  .FlagO .Flagl ,Flao2,F laqJ,F laj4 


I! 


m 

rot  >35  TO  SI 

C(I)«I 

READ  «iAkU>,l»(t>,Fi»(t> 
LW^KET'IHd- 


•FOR  All  Nil)  PRESS  Kll  FOR  TRACE  VARlAitEfi  Tail):  IP 


271 


?!! 


il 


PAUSE 

DEF  FNOd(A,B.C|*(-B-tSOR(r2-IIAICI)/2/<. 

DEF  FNStl«ufT]ifTii2tTk  I  Sm  gitum  Tl  atint 

KF  FMSi-Tii*C»+Oh*H2»*Kl+0»Hnc**Hco<Cft2i+2«(Ci2*Oi!+TJi2> 

Sm  qaaint  a  at  mi 


SI 


EF  FNUo4lia-FN3i 
FRtipi 


?7 


% 

J7I 

JHQ 

370 

410 

411 
421 
431 
440 
451 
460 
470 
ABO 

4S 

sit 

520 

530 

540 
550 
540 
571 
5B0 
570 
to  0 


DEF  FN j»in*Til *201  12+301 13+T 14 ! 
C  Ml 

W 

p*t-tn 

jttri ' 

DlO-.l  I 


6m  rlniinud  0 
Sum  Tl  tiidot 


Default  tital  nagi  It  0,1  kq/cu  it. 
Default  final  tiloM  U  1.0  ci  n, 


Djfaul t  nuUo  It  cinutnnt  «i1mi 
initial  noninal  valtat 


Fractional  irrir  f«r  c tnnuu rqirtcu 
Fractional  chonga  in  varlabli  to  cinpoli  dorliatlii 
Dlt  nay  bn  churigad  to  or  fron  ,5  by  Uj  KB 


Fat  It 


INPUT  'Nano  if  cinpiind'^O 

0I5P  "Firnula.  «  a  tins  par  null  «f  C,  H,  N,  0  In  irdir' 
'aw  k7  fir  TNT'j 
JT  Ac,  Ah,  An  .An 


Aifrvi  m,  >niiiNn  in*« 

M*t-M(2l>-*i(22»-Tic-N(19>»4l<2R)-Tin-H<43)«J4<44)»Cr-m23>«Ci 

Ct»Nu*Flaol»Fla(|>fltm2«Flao3«FU54«Cflao,Tinflaj»Ticflaa‘ 


Cl»1 

0 


IF  Ac*Q  THEN  AflE-ll 
IF  Si-1  THEN  ®«IE-10 


hM  THEN  , 

B1SP  'Internal  antroy  if  firnatiin,  J/nila'j 
DWJ  'j»o  SHIFT  k7  (-023)  fir  TNT') 

Fi«JLll21lll5fAh01.IIStAnll4,llb7*Aia»tS,7V?l 


Flipu  and  cindinud  phaut  nut  fir  mu  fill 
Pruldi  fir  full  with  m  C  ir  ni  H 


Firmla  nan,  q/mli 


c,ntl  (piMi1  ^ 


BMLMl  "  "Wl,  C/H',C.n 


20 

130 

til 

640 

670 

JR 

710 

710 

720 


Vpi 


Pl'BfTERM1  Mf 
PRINT  'Tital  cine  «  *)Cpn 
HtKpn/UtC  nil 
Nf>Cpii-HU  I 
CmNC«0  I 
Anp*PR0UNtlAn,-3> 
Acp*flOUNI(Ac,-3) 
Ahp'PIOUNDIAN  -3) 
Aip-ffQliiltA»«,-3) 
jJ*W 


Tital  natt  if  titanlin 
Hilt  if  fill 

Cointer  fir  itoratlim  and  crltirlm  fir  'Appru' 


Int'i  iltACE  / .  / . igC ini  \ iff'11 1 ^ i 1 ' i P '  H',K*  N'.K*  Q* ,K2)TC/t1-  M0Z.3D/ 
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«n 

821 

830  Cincai 

840 


Tihpi 


HMt'IWXUOOOMf/Fm.ir-iO) 

ite’ ' 

rmiji  It;*  Min  if  Ti* 

Bfcanr.'To.w Vn“ 

Nc*AtlHll* 

KiHo1mW,96  I 

N|a»A|aOI1oJs»iii.9!ill 

tli'l  • 

TK'l«t  I 

Al»»‘N'  ! 

VF  AIOO'V*  UttN  Tt«p 
T'.'SOI  I 
GOTO  Ttnjil 

IF  Flqgi  THEN  P1SP  ,:iutV33 
IF  Flail  THEN  DIBP 
ir  Flag2  THEN  P1SP  *Tn*2112 
IF  Flag?  THEN  CISP  *Tm2M7  ‘i 
IF  Fhai  THEN  DIM*  “’h«21«  1 
IF  Tlcrlag  THEN  PISP  *14-1291 
IF  I inflag  THEN  DISC  *Ta<3220  •; 

INPUT  "Appriniwti  nmirattrs  **.1 
INPUT  *#i  yiu  mm  'Appni'TU/NP.N 

ir  upcimn'r* 

II 

ii«*  ■ 


Then  nu>c 
THEN  Nu»l 


itJI  Tcnpl 1 1 
HU  Calc i 
1140 

HS1 
1  l/l 

nu 

n« 

'.200 
1210 
1220 
12.it 

1240 

U“J 

Ml 

£1 
13DD 
1311 
1320 


D1SP 

WVll’i  '• 

IF  Doll  THEM  it 
OISP  14 
fiyED  4 

"HINT  *1  ■  ■iTi"  N 
I'ISP  *T  •  "it,1  r 
DISP  D< 

vmm 
if  us-t  no  i» 

T2»T1 

U2-U1 

llrl 

IF  TJOI  THEN  Ii«p£ 

IF  Alt 0*Y*  THEN  fug 


T  lau  * 

/P*IV 

SFJ(V. 


•till 


Or 


J'  LIH2) 
‘ipf*  JO 


T-3SOO-SOU*SGH<IHi> 

1331  GOTO  Two! 

1340  T<hp2;  IF  A89(I»IUOOO  THEN  Ek 
IF  At»»V  THEN  Info' 

PI5J  Tf  * i i *  I  priii  Incj  nQuuCi  tinyirctfirs 
DISP  *  Inuroilam  cwf- 
1380  PAUSE 


13S0 

IM 


praiin'  an  ‘iNii’mUi* 
A^g  in 


Incl vdti  Nliti  of  N  atm  frm  air 

>  • 

Difailt  «alii.  TIC  ant  graghiu  not  tilth  prmnt, 
Ronti  lntirnai  mirov 
F w  tihMPatiPi  IMorpiUtnn 
Default i  nantai  unpiratiri  wtiing 

For  aitiMtic  tinoirgtiri  n'ting 


"rinding  giu  Mltmg  oimti  if  ctndmid  ghaut 


Plfaul*-  *T*  fur  niu  cincintratnni  *N*  itharulio 
Mpproi"  it  grnirailv  dm-obli  byliw  3010  K,  bat  nit  ab ova 


Cmoifii  oailhlirtim  cinitann  at  T 

rrtrtuimi  caltuhm  noli  nunbn  and  'iildaal  Polio  U. 


K4‘.*: 


Ati  mot  win  Niv  <5 


1390 

not 

141.0 

1420 

1430 

1440 

14S0 

1461 

1470 

1480 

1491 

1S00 

1S1I 


Intpi 

i 

En 


T  T7-U2d(T2"Tl);(U2*b) )i 
CulO  Calc 


Hmih«N*»5 
roo  1*1  TO  44 
H'«Nt*N(I)l 
HmnaNipniCUltNU)  ' 
NEX1  I 
Conc3;Ng«N»un  ' 

HaM-(MnHm.i692>*U00' 


Ha»g-Hais-E9 ,  n*Tip  -47 ,?»Hi t-63  9*T  i  i  -147  ,B»Ti2 
Hasg»ltng-223.7»U3-7".9*Ti4-41  9UT;n-12,0U*Cr 


'.moor  murpolatno 


Total  ailos  of  or oducov 
Hiltc  of  go; 

Hah?  c(  gao 
Ton'  nan  in  can; 


Nan  o r  gas 


‘an 
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1520  K„g*H«fl/Wq  I 

$11  WSfUffll'ldciniT 

1551  iMcmi _ MACE  /*T  • 


A*ga  Hilar  Mtt  •  f  gawa 
Aigt  Hilar  nh  of  all  pridaels 


m 

1581 

1S90 

m 

162 

163 


pif§T  *Aw  Hilar  nats<gaM»>  *  *)H  gi*  g/Ml*; 

PRINT  BPAT2)j*Am  Hilar  HwaltitalV-  *’Nj}“  g/«ol*i 

1 1>^**rpr*ftsurt  •  *IMI*  bar«*)LlNl2) 


VALUE  Id  BARS 


II  m 


.  it 

mi 


Z*)SPA(4)|*Part.  pm.*; 


1671 

m 

171  ' 

7 


II  A  V  tlVffWII 

1721  ihcint 


FOR  1*1  TO  34 

SWIAWBSUe-h 

PRINT  US' . 

nr.  m3 1 
mcSiNEXr  I 


Final  iitpit 


■i  inm  wupufTi 


I 


0 


U..B  8Al5X,i.3tt|5X,I.3DEl‘niU«‘ 
FOR  1*19  TO  23 
if  N(i)*t  Then  Cmc6 


C.«C6,  Mf'*****™*" 
1771  FOR  I«3S  TO  44 

17BI  IF  N(IM  THEN  Cane’/ 

liill  Cmc7i 

1811  Cupli.CMUl  jjhl 


1820  '""  PRINT  uSlNC  luciNu.Hf  iNl.AcpiAhPiAnp.Aap.C.H 
i$jl  Ittefi  MWT  ffi  ’P *l,r!4sfll!‘  Ourprmirt  ■  *,38, 3D,*  ban*, 


l! , 


lRIll 


PI  INI  LINI4) 

EP 

T  TOU  NISH  TO  RUN  ANOTHER  CONCENTRATION  <Y/N)?\M 

F  Ci«“Y“  THEN  Cant 
•Ana 


Veluna  ■  *,3D.3D,*  ca  h* 


INPUT 


•thir  full  (V/HI'jCrt 
jjjjj  gtl/fU|*i4^EN|YOFT|!^TI^ 

JJSJ  WUj^ljB.rROH  LCT(Kp>,PARAHETFIl|  . . ,  WWMfltMWiM Wl 


190 


1970 

1*?H0 


»  1*1  TO  SO 


lWO  LglmAklI.i)iAk(I.Z)/(TiAk(I,3))iAt(1.4)tT 

56*5  t.3k*LakiAk(i.6)0UT(U  rt>  I 

2010  lF  <1136)  AND  (1(41)  THEN  Lqi»l 


m 

m 

mi 


LqHgk/2 
N  Lgk*MNIlgk,V¥) 


m 


M-SP 

NEXT  I 
KticHXdV)  I 

|F  T|j®yi.  THEN  XUc-1120 
"  — THEN  XtH((22) 


Conmrilm  Tran  Ip  to  Kn 

Um  1>UR(K )  fir  U203  and  Ti305  lo  avoid  owfliu 


IF  T( 
K1*K<35 


. um 


THEN  Xl*k (36) 


2141 

Ml 

2170 

2181 

£81 

2211 

2220 


IF  T<2112  THEN  N2*X(38) 

X3«X(39)I 

IF  T(2I47  THEN  K3*K(40! 
i4*K(4i)l 

IF  T(2  43  THEN  K4»K<42) 

Xnn«X(43)| 

IF  1(3220  THEN  Xtln«4l(44> 
XI43»(K3/R4)*4/Rtl7K4 
M32*K2/K  t  IKK2/K3RK2I  ( X2/K3 ) )  *3 
K»31*Klt(Kl/Kti)»(Ki/Xtl>»(Kl/X3>ttl/K3 


Fir  TIC 

Fir  cimdantid  natal 

Fir  TiO 

Fir  T1203 

Fir  Ti30S 

Fir  Ti 02 

Fir  TIN 


224i 


I 


RETURN 


Thin  lait  3  mid  for  Lo*  ealan. 
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22$o  i  mmummtummmumnmmuutnmmtuummmmmt 
a  i  tmmtmMmnmhmunmmtmttMttttttmimmttmmm* 


II 

2301 

2310 

SI 

M 

8» 

23B) 

2390 

2400 

2410 

242) 

2430 

2440 

2450 

2460 

2470 

2481 

2490 

250) 

2510 

252) 

2530 

2540 

2550 

25611 

2570 

2580 

2590 

2600 

2610 

2629 

2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
27S0 
2760 
2770 
2780 
2791 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2*00 


C«=0  ! 

IF  Nn  THEN  Nut 
Ncc*Nc 
tocsin 


C«»nts  iterations 


FUR  I-i  TO  48 
N(I1«0 
NEXT  I 

Ar»N(6>«.4036 

fl119II5fla8'*f[1®92’,rls93=floo4*Tlnflag*l'lCT'lanaL«sfl69e0 
ZnfiRiNn/e)  !  Applies  when  Tinflag“0 

IF  NoalNc  THEN  Apix 


Co*N»o  I 
Ticflugal 

IF  NtilNcc-Hea  THEN  Aoptt 

Cflcg*i  1 

Act"! 

X=Ca/Acc/K<7) 

Ti«N(i)*i/Acc/!UlC 
GOTO  Apqx 

AppmFlagOci  I 
Ti*l/KU 
Acc*l/Tl/I(tic 
x*ro/x/oi(7) 

GOTO  A^px 

ApsxiCe*Ncc  I 

IF  Nti<Nw-Ncc  THEN  Appl 

F lag l*i  I 

Tinfiaoxl 

E<U«Ntl-He«-Ncc-Nn  I 
IF  FxtilO  THEN  Npp2 
IF  X rin/X t tK2/Klil(2/K  1)1  THEN  Appia 
Z3SBR ( -Ex 1 1/2)  ! 

Ti*NllM/Ktin/Z 

I 

App3:X=l/U/Tl 
GOTO  Appx 

I 

App2:rlag0«l  i 

THRlNIExtl.l/IU’ 

Z*l/Knn/Ti 
GOTO  App3 

ApplilF  Nti<2*(Noo-N«>/3  THEN  Aop4 
X=K1»KVK2/K2  ! 

Flanl«flag2»I 

Ti«f/Xl/X 

IF  KtIn/KlM2/Ki*K2/Ki<l  THEN  Aupx 
Tinflag1! 

Z*l/Ktln/Ti  I 

IF  Nu-2t(Noa-Nci:)/3>Nn  THEN  Appx 
If  64LOC(K3HMUnnl>10»LQG<F2)  THEN  Aoo4a 

Appla : 2=30R <  <  Nn-Nt i *2»(Noa-Ncc 1/31/2) ! 

Ti*l/Ktin/Z 

X*(l/K2/Ti)*(2/3> 

GOTO  Appx 

i 


No  oxides 


Graphite  present,  no  natal. 


Ne'al  present,  no  graphite 


Natal  oadis  present 

Excess  li;  Ti'TiO 

The  excess  of  Tj(c)  over  N 

N  is  in  excess:  TiN  +  T.O 


Titc)  i;  in  excess  oner  N 


TiO  Ti2U3  region 


TiN  *  2  oxides 


TiN  *  T i203 


32 
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m»w  trom « «» « 


i  PrtflgilF  TiAflM  THEN  PRINT 

■  <r  v  i  nun  in 


IF  Ticttoo  THEN  PNXKT 


TIN*  | 
TiC’j 


IF  Cflafi  Ad  (Ac >iE-iO)  THEN  PNIMT  *  6r*i 
IF  LixflM  M.LflM  THEN  PRINT  *  L«*j,(l;Lfj‘ 

if  Fisfi  Then  print  ■  n*i‘; 

It  Cl. It  TUtU  ItlUT  •  Tin«. 


...  THEN _ 

IF  Flag!  THEN  PRINT 
IF  Fla|2  THEN  PRINT 
IF  Flog3  THEN  Ft‘“ 

'TOT 


)■) 


INT 

PRINT 


Net* 

no* 

W\ 

pimn^"  . .  Tl0Z  1 

RETURN 

MtutmMwmnuttnmwutiMttmwuuMmttMMM 

DISPLAY  FLAG  STATUS 


ipflgi  lisp  TABU)  ]  SliplQp  flag  itQtus 


3721  Dtpflgi  DISP  TABU)  i 

3731  DISP  “NEKTON'S  CALC.,  •lSPA<S>t 

T740  IF  Tlnflag  THEN  RISP  *  TiR'i 


aa  THEN  DISP  *  TIC* 

IF  Cflao  AND  <Ac)iE-m  THEN  DISP  * 
IF  LikFIm  OR  Lttag  THEN  DISP  *  Lex1 
IF  Flag#  THEN  DISP  *  Hit*! 

IF  Flagl  THEN  DIBP 


Gr*j 


(‘iLfi*)*; 


IF  Flo#  THEN  OISP 
IF  Flaq3  THEN  DISP 


THEN  DISP 


TiO’i 

T12C3*; 

Ti30S*j 

Tia2"i 


MSI  Ncwti  I  A  NEWTON'S  HETHOO  FOR  UP  TO  FIVE  SIMULTANEOUS  H0IHIME«  EONS 

mi  I  wiHttimmmHiPtitmmtt»ttiHttittM«iniiiHtuuunuu»u 

3jS70  REDIN  X»(II),Y«(W>, Will), WUli), Dill, 

3080  COSUB  DtpFlg 

1891  NeutOiCv*!  !  Oomtir  f»r  Neuter's  Mined  itiratiens 

Huct*Nii(t+l  I  Cienti  entries  inti  Neut 


(0  IkutJi 


3920 
1931 
39(8 
3951 
3961 
3970 
3980 
3991 
4008 
4110 
4020 
4031 
4040 
40S0 
4060 
4170 
4IB0 
4090 
4i  0 1 
4110 
4421 
4130 
4141 
41511 
4160 


bct+1 1 
HAT  xt»: 


Xn 


PRINTER  IS  7.1 
Cs*Cs*l 
FOR  1*1  TO  li 
K«1 

COSUB  Fx 
YnU)*Fx 
FOR  J*1  TO  II 
X*J 

Xt(J)«(l-Dlt)«AtJ)l 
CO  SUB  Fx 

D!l,J)»-(Fx-Yn(l))/Xn(J)/Dlt  I 
Xt(i>-Xn(J) 

NEXT  J 
NEXT  1 

IF  DETjDJOO  THEN  Ntwtl 


Dlt>.i  by  default:  nay  bt  resst  te  ,5  by  K6 


Partial  »f  Tn(i)  urt  Xn(j) 


PRINT  “MATRIX  SINGULAR* 
[  SINCUUw* 


DISP  “MATRIX  1 
PRINTER  IS  I 
PRINT  LIH«1),*D  Halrix;* 
HAT  PRINT  D 
PRINT  UNU) 

Sun 


1 


4171  HtutliHAT  Iti*WV(D> 
4180  HAT.  Dili*InlYn 


Aire  HAT  Xii* 

4200  N«ut4 1 FOR  1*1  TO  Ii 
4210  IF  Xn(I>)*0  THEN  Neut2 

4220  Xn(l)*()bi<lHDch(I/)/2! 

NEXT  I 


4230  Neut2: 


If  <*  0,  use  1/2  previous  value 


4240 
4250 
4260 
4270 
4281 
4291 
4300 
4311 

$30  Retn: 


FOR  1*1  TO  11 

IF  ABS<Yn(I)))EbS  THEN  N«wt3 
COSIIB  Ft 
NEXT  I 

PRINT  Cut"  ITERATIONS  FOR  “;Ii}“  VARIABLES! 

GQSUJ  Prtflg 

PRINTER  IS  7.1 

IF  Chart*}  THEN  RETURN 

GOSUB  Sun  ! 

RETURN  I 


Dugnostics  given  when  Chcck-0  (given  by  key  XI)) 
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4340  i  mmmmmmutmmutumuutunimmunnmmnmutm 

4371  i  Sets  flaps  far  condensed  theses, 


Tinl*Tinfiag  ! 

IK.fi  Of 
L**C-Lo*f laq 
FOO'FlapO 

4459  $041$ 

F3l*Flag3 
F4041og4 
Tlcl-Ticfiaq 
CO-CflM 

CflegvHlNil.INT<Acc)>*<Cr>*l> 

IF  Gr<»  THEN  Rr-I 

IF  AcME-U  THEN  Cflao-i  I  ,  . 

TKna9*(Ti»tecll(Uc>l-iE-5)IITic)*»)*(4£)lE-it) 
IF  Tlcfl  THEN  Tlc»l 
lf«NOT  (Tlcflaa  ^44D  Cflaq) ' 

4570  Tlnflan*imzHiin>l-l£“S)4<Tin>«|t 

4580  IF  Tln((  THEM  Tin-4 

4500  Flaql*(Ti»Kti>HE~S>*llWt>*l)i 

4401  IF  H«t<0  THEN  Hti-4 

m  re<?m-irs,*,TU>‘D> 

4430  Flog2«(T10XtSOI(X>>XI;)l-lE-5>  4(112  >*9 ' 

4440  IF  T12U  THEN  112*0 

m  fmmyp™1™-™'1'*™ 

4470  Flao4*(Tl4X4X4K4)l-i£-5)4(Ti4)*0) 

4480  IF  T14(0  THEN  T14-0 

4400  COSUB  Act 

4700  LofIag*4T)2143)  AND  (U*»4>  AND  (AtsutOl-Ens) 

4710  I 

4720  i 

4730  IF  Lox(0  THEN  Llx-0 

4740  IF  l.ixflag  THEN  Flaa0'=Floqi:*Flag2«FUa3*Flaa4*0 

4750  Lflao-Loxf lag  M40  Lf 

4'40  IF  Lflag  THEN  LoxflagMIi 

4770  I 

4780  Tc«R=Nt  x-FMSt i-T m  ! 

47V  0  OconMIoa-FNSo  ! 

4800  IF  NOT  (Cflaq  AND  Flaqi)  THEN  Flaga 

4910  Cflag»(Ncc)HlN(Tcon ,Qcin)/3)  I 

4820  FToalM-Cflag 

4830  FliqEMicflOdM 

4940  KlagailF  N01  (Tlcflag  AND  Flag3)  THEN  Flagli 

4650  Flqq2*CflagM  ! 

4840  Flua3-(Tc(»n)i0*Mcc) 

4B70  Tic? lag-i -F laq3 

4880  Flaab : JF  NOT  (TUflag  AND  Flag4>  THEN  Flaga 

4890  Flaq2*Cflaa*l  I 

4900  Tlcflag*(2iOcen/3(Tcin) 


ForMtiin  constants  post  be  exceeded 
and  apovnt  <f  phase  nit  negative 

1  a  save  tld  mint 


Far  fails  with  no  carbon 


IfM  unless  TiC  and  graphiti  both  present 


For  Ld  ;  Sun  of  activities  if  liuvld  oxides  +  petal  >*1 
anil  T  >  mlt,  pt,  cf  T102 


Lflao  for  Lox  exesat  Lixflag  us«d  when 
both  TIC  and  graphite  also  oruent, 


4911  Flaq4»i-Ticf lag 

4920  flagciIF  (Tkflag*l)  UR  <Flag0*l>  OX  (Cflag»t>  OX  (AcME-10)  THEN  Iiset 
4930  FlggO=(Nti-FN3ti)Hcc>‘  Ni*h  TiC,  timer  Net 

4940  CflagM-flagO 


both  TIC  and  graphite  als 
Condensed  Ti,  except  for  TiN 
Condensed  O 

Or  f  Til  *>  TiC  t  Ti2 


TiC  ♦  Ti3 


TiC  *  Ti4  »>  C  *  TiC 


4940  CflagM-flagO 
4VS0  ! 

4944  Uset;Iiflag»Flag0tFlagl«-Flag2tFlag3tFlag4tLflau 

4^80  '  IF  Iiflag)2  THEN  Flagl 

499  0  1 1 =4- Iif  I.qg-T  inf  iao-T  icf  iag-Cf  lag ! 

SIJ0  ir  11*0  THEN  RETUXN 


or  Gr,  nit  both 


I iflOQ  can  be  0,  1,  or  2 

Iifloo  =  1  for  lflag  but  0  for  Loxflag 

Nupber  of  variables  in  •Newt’ 


RED1H  Xn(h) 

1  Following  sets  values  of  variables  in  ’Newt* 

IF  NOT  T inf lag  THEN  X»(I1W 

IF  NOT  Tinflaq  AND  NOT  Cflaq  AND  NOT  Ticflaq  THEN  Xntli-l>*Acc 
IF  Tinflog  AND  HOT  Cflaa  AND  NOT  Tlcflag  THEN  Xn(li)*Acc 
IF  Iiflag(2  THEN  XndkX 

IF  (IiflagM)  AND  Tinflog  AND  NOT  FlugG  AND  NOT  Loxflao  AND  NOT  Lflag  THEN  XnUiMi 

IF  (iiflog=0)  AND  NOT  (Cflaa  AND  Ticflagt  THEN  Xn!2Mi 

RETUXN 

The  variables  are  X  ■ if  no  pore  than  one  oxide  or  petal), 

Ti  (if  no  oxide  or  petal)  and  Z  (if  TiN  is  absent). 

Reception;  When  TiN  and  one  solid  oxide,  but  no  petal  or  liquid  oxide 
cresent  then  Ti  (and  not  X'  is  firs*  variable. 
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5145  Fhgi!Tcon»Nti-FNSti-Tin-Tic  1 

5149  IF  Flogi  THEK  F13a 

5171  IF  NOT  Flag)  AND  NOT  Flan  4  THEM  F13b 

5180  IF  NOT  Flag3  MO  NOT  Fla*4  THEN  F13c 

5199  IF  OconHcon  THEN  Fla  ! 

5219  Flag  0-Flag  1*1 

5211  Flag2=Flaq3*Flag4=0 

5225  GOTO  I  it<l 

5235  Flat  Flag0>4 
COTO  F14a 


The  cast-  of  Iiflag>2,  which  is  disallowed 


The  case  if  Iiflag=5 


5244 
5251  F14i 
S260 

5274  FUd 
5280 
5295 
5349 

5310  F13b; 

5324 

5339 

5344 

5351 

5360  ri4ai 
5370 


IF  NOT  FlagO  THEN  F14e 


Fir  Iiflag  *4 
Flags  0,1,2 


Flaq2=IOcon)Tc»n! 

FlugO-NOT  Flog2 

COTO  litet  rl  ,  „  , 

i  Flags  l,i,3 

Flag9*Flan4=0 
Flag3*(Ocin)l.tOTcin) 

Flaji»N0T  Flag3 
COTO  list! 

!  Iiflag  =  4,  FlagO  -  0 

IF  Qcon(=l,S4Tcin  THEN  F13b 
Flagi=0 

t  Float  2,3,4 

Flaq4=(0cin)54TcMi/3) 

Flag2«N0T  Flag4 
GOTO  list? 

5431 !  mmmttmutummnmmmmmmtmmmttmmnmmmu 
5440  Fx : !  054000004  Fitting  functions  far  Hinton’s  nethod  calc.  04000000000000 
5459 1  tmmmuuuittuiumummmnMumtuumtmmummm 


bJ/O 

5389 


5399  F13a 
5449 
5419 
5429 


5469  1 
5479  ! 


i 


I 


5464 
5491 

5599 
5510 
5529 

5539 

5540 
5559 
5569 

5574  FootiliOSUB  Biff  ' 


Ft*  conoutts  errors  in  stoichiwetric  conditions  for 
exisitinq  ualots  of  variables, 

Finds  appropriate  values  of  master  oariables 
Calls  on  ‘Sotc"  'a  cortoute  nelcs  of  all  snecits 
)tt(  i)  is  used .  since  both  oriqinal  value  if  variable 
and  value  changed  bv  Bit  are  enpUyed. 

IF  NUT  Unflag  THEN  Z»Xt(Ii) 

IP  NOT  Tihflao  AND  NUT  Cflag  AND  NUT  Ticflag  THEN  Acc=Xt(I;-l) 

IF  Tinflag  AND  NOT  Cflag  AND  NOT  Ticflag  THEN  Acc=XTUi> 

ON  Iiflagti  GOTO  FxO,Fil,Fx2 


5580 
5599 

5609 

5610 
5629 

5639 

5640 
5651 
5661 
5679 
5689 

5691  Fid: 
5790 

5719 

5720 
5734 
5749 


IF  <I=1)  AND  (IifloglS)  THEN  F).*Do/Noa 
IF  (I«l)  AND  (lifiag*!)  AND  Tinflag  AND  NOT  Lflag  AND  NOT  flooO  THEN  f x^ti^N’i 
IF  Loxflag  AND  ( < 1=2)  OR  (1=1)  AND  Cflag  AND  Ticflaa)  THEN  ri=Atsun-l 
IF  U*2)  AND  NOT  Iiflag  AND  NOT  Lflag  THEN  F.t=Dti/Nti 


Exit  routine  for  “Fx" 

Kollauina  are  stoiChiiHetric  conditions  to  use  in  'Newt 


IF  NOT  Tinflag  AND  ( I=!i)  THEN  Fx=DeIn/Nn 
IF  NOT  Tinflag  AND  NUT  Cflag  AND  NOT  Ticflag  AND  (I=Ii-i>  THEN  F»=Dt)c/Nc 
IF  Tinflag  AND  NUT  Cflao  AND  NOT  Ticflag  AND  (I=[i)  THEN  Fx=Delc/Nc 


The  stoicluonetnc  tests  correspond  with  the  variables, 
Returns  to  "Newt'  (or  to  "Nw*  for  li=0) 


2) 


For  ltflw=0 


RETURN  ! 

i 

X=Xt ( 1 ) ! 

IF  NOT  (Cflag  AND  Ticflag)  THEN  li*N(l)=Xt( 

IF  Cflag  AND  Ticflag  THEN  Ti=N(i)=l/Ktic 
IF  NOT  Lnxflog  THEN  Foot 
G0SU8  Act 
GOTO  Foot 

5759 1  tmmmmtmmnmtmmmmmmmtumittummtttiitu 

5769  Ad:!  OOUOOOOOO  Gives  activities  in  liquid  oxide  1196  0  909090909*000999909 

ttmmummtmtiutmmtmmmmmmnmttmtmntumtt 

Atl=KtiOTi!  Activity  of  metal  in  liauid  oxide 

Following  arc  activities  of  oxides  in  liavid  onde 

AtlnRlOTlOX 
At2=TiOTiOX9XOXOK20K2 
At3=TiOTiOTilXOXOXOXIXOX30X30iLflag<2' 

At4=X49TiOXOX 
Atsv«»At9t  At  i*A  t2*-At3+At  4 
RETURN 

mttuuhmtmtumutmnmmmmmmtmttutmtummtt 


5779  ! 

5789 

5799  i 

5899 

5811 

5820 

5830 

5840 

58SI 

5860  ! 
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5870  Fil: • 
5889 
5891 
S90I 
S919 
5921 
5930 
5940 
5950 
S%0 
5971 
5980 
5990 
619  8 
6811 

m 

6140 

6959 

6069 

6870  fxic. 
6089  ! 

6690 

6108 

6119  Fxlbi 

6120 
6130 

6149 

6150 

6169 

6170 
6139 
6190 
6200 
6210 

6221  Fxla: 

ill 

6250 

6269 

6270 
6230 

6290  r x2: 
6300 
6310 
6320 

6330  0x2(1: 
6340 
6350 
6360 

6370  fx2b: 
6380 
6390 
6400 

6419  F  x2c : 

m 

6440 

6458  Fx2d: 
6460 
6470 
6480 

6490  F«2t: 
6510 
6510 
6520 


IF  T  inf  log  AND  HOT  f  lagi  AND  NOT  lflag  THEN  File 
X=Xt(l> 

IF  NOT  Lflag  THEN  Fxib 


Fir  Iiflog  *1 


F or  Lflag^l,  solves  9  cubic  t9  find  Ti  =  f(X) 
Must  satisfy  sun  »f  activities  =  1 


Cotff.  of  cubic  (Unit  cutff.  For  cubic  term. ) 

Reduced  cubic  lacks  quadratic. 

Discrininant  of  the  cubic 


Avoids  cube  roots  of  negative  no. 

The  onlv  real  rut  af  the  reduced  eqn. 


A3=-t/K3/X3/X/X/X/X/X 
A2=-(Kli*Ki»X4k48XIX>*A2 
Ai=-X»X*XtX2«2*A3  ! 

Ar«<3*A2-A18Al>/3 
Br=(2Mittl«Al-98Al*A2)/274A3  ! 

Dt=Brttr,'4+ArMr»Ar/2?  ! 

IF  Dt(=(l  THEN  Fxtc 
Ri=-Br72tSW(k) 

Rl=SGN(RM(SGN<Ri)Ml)Ml/3i  ( 

R2*-Br/2-SQR(fc) 

82=SCN(82)t(SCN18210R2)Mi/3) 

Xr*Ri+R2  • 

Ti*N(i)“Xr-Al/3 
IF  Ti)0  THEN  Feut 

IF  Lflag  THEN  Ti=N(l)=FN9dIX8XtX«X2«2,Kti+XWX+K4tX*X.-l) 

In  case  cubic  bad  negative  r oo 
Lf=2 

GOTO  Fout 

i 

iF  Flaqi  THEN  Ti=N(ll=i/<ti 
IF  Flag!  THEN  Ti*N(l)=l/Kl/X 
IF  Flag2  THEN  Ti=90(  1  )=1/X2/X/SRR(X) 

IF  Flaq3  THEN  Ti*N(l)=(i/X3/X/X/$»<X>)‘<2/3) 

IF  Flag4  THEN  Ti=N(t)=l/K4/X/X 
IF  NOT  (Cflag  AND  Ticflaq)  THEN  Fout 
I=Xt(i)l 
Ti=N(l)=l/Xtic 
X=ti/X2/Ti)*(2/3) 

GOTO  Feut 
ri=N(l)=Xt(l)  I 
IF  Flaal  THEN  X=l/Xl/Ti 
IF  Flaq2  THEN  X=(VK2/TD*(2/3) 

IF  Flaa3  THEN  X=(l/K3/Ti/S08<Ti)>\4 
IF  Flag4  THEN  X=S0R(1/K4/Ti) 

COTO  Four 


Ti=N(i)=l/Kti  ! 

X=l/Ti/Xi 
COTO  Fout 

IF  NOT  (Flaqi  AND  Flag! 

X=K1A28K1A2  ! 

Ti=>N(l'*lAi/X 
GOTO  Fout 

IF  NOT  !Flaq2  AND  Flog 
X=<K2/K3)*49K2tK2  1 

Ti=NU)=i/X2/X/SBRCX> 

GOTO  Faut 

IF  NOT  (Flag 3  AND  Flag 
X=X3/X4A4A48X3 
Ti=N!l)=l/X4/X/X 
GOTO  Fout 

IF  HOT  (Fhg2  AND  Flag 
X=K2A4IK2/it4 
Ti=N(0=lA4/X/X 
COTO  Fout 

IF  Cflaa  AND  Flagi  AND  NOT  Ticflaq  THEN  Alert2 

Ti*iAtlc 

X=l/Kl/Ti 

GOTO  Fout  ! 


Cose  of  lit  ♦  graphite  +  Ti203 


Fen  TiN  +  re  oxide 


THEN  Fx2a 

Ti  -  Ti0  region 

THEN  Fx2b 

U0  -  Ti20C  region 

THEN  Fx2c 

Ti203  -  T 1305  region 

THEN  Fx2d 

THEN  Fx2e 

NWCTP6544 


6330 

m 

6561 

6590 

6530 

6S96 

6600 

6611 

6620 

6631 

6640 

6650 

6669 

6671 

668') 

6690 

6700 

6711 

6720 

6730 

6740 

6750 

6760 

6770 

6781 

6790 

6801 

4810 

6820 

6830 

6840 

685C 

6860 

6870 

6680 

6890 

6900 

&9lO 

6980 

6r'30 

6940 

4950 

6960 

l9r,0 

6980 

6990 

7000 

■’010 

7020 

7030 

"'MO 

7050 

7060 

7070 

7080 

7090 

7100 

7110 

7120 

7130 

7140 


IF  Tinfht  THEN  Z=l/Ti/Ktm 
IK  Ticking  THEN  Acc=l/Ti/Ktic 
50SUR  Spec 

Tin=(Nnc-No-2<N8)*Tlnfl4g 
Tjc-(N:c-Ci-Ci2)ITicfiflfi' 

i 

Hi44)*lint(T(I226 
N!43'«Tih-N<44> 

N(20)*Tic*(T<3290i 

Ml?)«Tit-N(») 

COSUB  Oxides  !  GIVES  OVIKS  AND  ItSTALS 
Dn*Nti-FNSti-Tin-TiC-t1tt-TNJ‘ti* 

D»=N#i-rMS«-Tli-3»Ti2-5*T*3-2*ti4 
De )  n«Nnc*No-2IN2*T  ir. 

IF  Cflao  AND  Ticfkg  THEN  T*c*im>Tjc 

Deic=Ncc-Cc-Ci2-Tic 

Gr*NCJ)*l>»U»Cn*««c)lE’ll* 

RETU»« 

i 

!  tmtmmnmmmnmuMuummuumtmuttmmmttmmt 

!  41*444144414*414  Connotes  ftclcs  of  nil  ooseoi1*  soecios  tltUluillOlit* 

mmtmttntuuitnmtummmmuiniutununtitntuuttmt 


*iic  gad  Ncc  or«  coe.-ected  nui  if  n  ir  c 
sis«  r,  gmn  ii>  'Seec’ 


Spec 


Nit: 


Co.-: 


•V  !'««;«  with  no  co^Oon 


1  leapt:  V  tj  1.  Ac: 

02«N<18>*X'2, 

0-N(17)*K(17)iX 
Tlo*M2>*K(2)tTitX 
Ti»2*N(3'sK(3>*T.*X*X 

N6-N(16'=3‘2 
No*N'15)«K<15>»X*Z 
tloi*«NU2>»K«2'IHg**X 
IF  Cf:ttc*l  then  Acc*i 
IF  Ac=lE-10  THEN  Acc=lE-30  ! 

Co«N'.7)=k(7)*w.c4> 

C«7=N(0'ek(8'«Acc*X*X 

Hhn(9i+kaO)»X+Acc*(X(2?)«4X(30)i24>nt3r«xi+»  Zi/fAcctAc: 

‘.6<fflC.il't  ti  (tOitOMi-OMH  StecUt 

DH-'.H'  12tW'"(i32)MctlX*K(33i4Acc*Acc  '  CuYYigHi*  jf  ijjhv'frgnen  spccitc 

v*FNfld(2»Ph.Mh,-Nh> 

h»N'9'=K9)*Y 

Oh=Rf  1C»*KMD.’*X*Y 

82o=  N( 1£'«K { 12 )iX4YTY 

Hcn=N<29>*K(29>*Acc4Y*Z 

Hna=e'30)=K(30)*Acc4xm2 

Hco=N(3i)«Ki31)ttcc*X»Y 

Ch2o-N(32)=K(32)»AcciX)fY*r 

'I2h2=N'.33)*X(33)$Acr.1AcctYIY 

Tio=N(4)=CnH*N<24)=S«R(K(4)»K(24)|Tj»tec»Z' 

CnsN(2S'=>K(2S)4Acc4Z 

C2h=N(26)=XI2A>IAcc»Acr|f 

C2n=N(27)sK(27)*Acc«4cc4Z 

C3=N(34)=K<34)4Acc4Acc4Acc 

Ncc*Nc’hcn-Hnc4-8co-Ch2'Ki'-Cn«-2t(C2h2+C2h>C2i>  ■  -340 

Crf'ec’id  *oles  of  C  n’-pi! 

Hic=8n-H':n-Hnc!i-|.n-E2n-Cne  1  .sr-ectfd  iclos  of  N  .roes 

IF  Nnc<0  THEN  Nnc=lE-3( 

SETURN  | 
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7iso  i  mMtmtmtmttmmutmmtmmnttnttnmuumtmmntnt 

7169  !  tumtttu  Hales  of  als  and  oxides  liJMlMOiimWWnattmi. 
7171  !  tittittttitumiumuimuummtmmuuniimtfmmmimiB 

71B0  Oxides!  FOR  In=K  TO  22 


7191 
7200 
7211 
7220 
7230 
7240 
7250 
7260 
7270 
7280 
’29  0 
7300 
7310 
7320 
7330 
7349 
’350 
7360 

7370  Oxb  l 

7380 

7390 

7400 

7410 

7420  Oxc: 
7430 

7449 

7450 
7*6! 

7470  flxd: 

7*80 

7490 

'500 

7510 

’520  Oxe: 
7530 
7540 
7550 

7560  0x1 i 

7SV0 

7580 

7690 

7600 

7610 

7620 

7630 

7600  Qxtin: 
7650 
’66! 

76’3 

7630 

7693 

7700 

7710  Oxlex 

7720 

7730 

7740 

7750 

7760 

7770 

7780 

7796 

7800 

7310 

7820  0x2: 

7830 
7840 
78S0 
7866 
7870 

7831 
7890 
7909 
791 C 
7920 


Fallowing  far  tut  oxides  or  nets!  +  one  oxi' 


H(in>=0 
NEXT  In 

FOR  ln=3S  TO.  44 
N(In)=8 
NEXT  In 

Het-Til=Ti2=Ti3=Ti4=3 

1  Preced ing  orevents  carrvover  of  old  values 

IF.  NOT  1  if  lag  AND  NOT  Uxflag  THEN  RETURN 

Tcon=Nti-FNSfi-Tin-Tic  !  Condensed  Ti  +  oxides 

Ocon=Nia-FNSi  !  Condensed  D 

IF  Loxflag  OR  Lflag  THEN  Oxlox 
IF  liflag-i  THEN  Oxl 

!  Following  for  tut  oxides  nr  ttetsl  +  one  oxi' 

IF  NOT  (Flagt  AND  Flag!)  THEN  0x9 

Tii=Occn 

Net=Tcin-Til 

CPTU  0x2 

i 

IF  NOT  (Flagl  AND  flag2)  THEN  Oxc 

Ti2=Qc«n-Tc«n 

Tii=Tcen-2«Ti2 

GOTO  0x2 

i 

IF  NOT  (Flaa2  AND  Flag3)  THEN  Oxd 
Ti3=2*Ocon-3ITcen 
Ti2=Ucan-3*Ti3>/2 
GOTO  0x2 

i 

IF  NOT  (Flaq3  AND  Flao4)  THEN  Oxe 
Ti3=2*Tcon-ucon 
T,4=Tcon-3*Ti3 
GOTO  Qx2 

l 

IF  NOT  Flag2  AND  Fiag4  THEN  Alert3 
Ti2=2tTcin-Ocon 
Ti4=Tc«n-2*U2 
GOTO  0x2 

IF  Tinflca  THEN  axxin 

IF  FlaqO  THEN  Net=Tcon  !  cor  hflac  =1,  no  liquid  oxide 

IF  FXccl  THEN  Tii=lcon 

IF  Flqg2  THEN  Ti2=Tcon/2 

IF  FiagJ  THEN  Ti3=Tcon'3 

IF  Flag4  THEN  Tin=Tcon 

GOTO  0x2 

i 

IF  UaaO  THEN  Het=Tcon-Nnct-No+?tN2  1  RovUne  for  TiN  *  one  oude.  This  hre  is 

IF  n'.noi  THEN  Til=Ocon 

IF  c lag2  THEN  Ti2=0c«n/3 

IF  Fiag3  THEN  TU-Odin/i 

IF  Flaa4  THEN  Ti4=0con/2 

GOTO  0x2 

:  S«hox-Tcon/U+Ar2e24At3)!  liauid  oxide  oresen’ 

IF  Cflag  AND  Ticflag  THEN  Si'noxsOcon/<i»2*At2*4*At3*At4i 

1  For  this  case  use  O  -oiber  than  Ti  nalance 

GtlSiJO  Act 

Nc'rfttCISssex 

Til=AtUSu««x 

Ti2=Al2*S»itox 

Ti3=At3lSu«ox 

Ti4=At4lS*nox 

Lox -He t *-T il-Ti2+Ti3+T i4 •  Holes  oJ'  licuid  o>ide 

N(22)=Het»(T(l533' 

N(2ii=Net-N(22> 

N!36)-Ut*(T(2623> 

Nt35>=Til-N(36' 

H(38)=ti2»(K2112) 

N(37)=Ti2-Nt33) 

N(4Cl=Ti3*(T{2647l 

N(39'=Ti3-N(4|) 

N(42t-Ti4*(!(2143- 

N(4ii=Ti4-N!42! 

’ETURN  i 


cor  hflag  =1,  no  liquid  oxide 


Routine  for  TiN  *  one  oiide.  Th;;  ijre  i5  ror  oqi  in  N. 
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7930 

m 

7960 
7971 
7980 
7990 
8001 
8010 
8020 
8030 
8041 
8150 
8060 
8070 
9080 
8090 
8100 
8110 
8121 
8130 
8140 
8150 
9168 
8170 
8181 
8190 
8200 
9210 
8220 
8230 
8240 
8250 
8280 
8270 
3280 
9290 
6300 
8310 
8320 
8330 
8343 
8360 
8360 
3370 
8380 
8390 
8400 
8410 
8420 
8436 
8440 
845  0 
3480 
9478 
8480 
8490 
8500 
8511 
8520 
8530 
0540 
8558 
8560 
3570 


*  mutmtmttMtmmnmmmmtmttttmtumnmmttmnn 

!  mtuutm  Energy;  Comoles  Delta  l!  tmumtMttWlUtMm* 

&iergyiDu=-ui 

FOR  1=1  TO  44 

U(I)*lotU)*8#(I.2)»T-»B»(I,3)mMua,5)/Tt8u(I.4)«1.0G(ri 

dolor  internal  energy  nf  each  siiecies 


9l*Dv6U(I)tN(I) ! 

NEXT  1 

If  Du<=0  THEN  RETURN 
IF  FlauO  AND  (T=1933>  THEN  EnO 
If  Ticflag  AND  <1=3296)  THEN  Entic 
IF  Final  AND  !T=2023)  THEN  Enl 
If  Flag2  AND  <T=2112)  THEN  En2 
IF  Flag3  AND  (T=2I47>  THEN  En3 
IF  FlagT  AND  (T=2143)  THEN  EM 
IF  Tinflag  AND  032201  THEN  Enun 
RETURN  ! 


Sms  total  internal  energy 


Entic;  N(20)=Du/71145  ' 

N(19)=NAX(«,Tic-N<20>)  ! 
IU'=Dut(N(l9)=0) 

IF  N(l?)  THEN  RETURH 

N(20)=Tic 

GOTO  Enbeep 

EnO:  M22)*Du/lB623  ! 

N(21)=HAX(I).tet-N(22))  ! 

Du=Du*(N(2i 1=0) 

IF  N(2i)  THEN  RETURN 

N(22)=het 

GOTO  Enbeep 

Emin:  N(44)=tu/64968  1 

N(43.'=NAX(0,Tin-N(44))  ! 

Du=Dul(N( 435*8) 

IF  N<43)  THEN  RETURN 

N!44)*Tin 

GOTO  Enbeep 

Eni:  N!36)=Du/54405  ! 

N(25)=HAX(0 ,Til-N(38))  ! 
Du=Du«N(35J=0> 

IF  N<35>  THEN  RETURN 

N(38)=Til 

GOTO  Enbeep 

En2  N(38)=Du/li0484  ! 

N(37)=tlAX(0  Ti2-Nl38) !  ' 
D«=Dol(N(37)=0) 

IF  N<37)  THEN  RETURN 

N(38)=Ti2 

GOTO  Enbeep 

En3:  N<40)=Du/i38105  i 

N(39)=HAX(0 ,Ti3-N(40) )  : Liquid  “1305 
Du=Dn*(N(39)=8) 

IF  N(39)  THEN  RETURN 

N(4fl)=Ti3 

GOTO  Enbeep 

En4:  N(42)=Du/85960  ! 

N(41)=NAX(0.Ti4-N!42'!  ! 
Du=Du*(N(41J=D) 

IF  N(41)  THEN  RETURN 
N(42)=Ti4 

i 

Enbeep;  SEEP 

BISP  "Select  a  lower  tenperature" 
RETURN  ! 


Return  frert  "fical"  to  rain  pna^an 
In  the  Pollening,  2  phases  ef  *he  sub; 
Solid  is  aivtn  by  dU  divided  by  energi 
Solid  TiC 
Lieyid  TiC 


Sell.!  Ti 
lionid  Ti 


Selid  Tin 
Uovid  TiN 


Solid  TiO 
Liatud  TiO 


Solid  Ti203 
Liavid  Ti203 


Solid  Ti305 


Solid  Ti02 
l lauxd  Ti02 


;Tance  are  present, 
■  of  fusion. 
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Kfct  !  utmmmmntumttummuummmtmtuttmtmmtmmtt 


mi 

841  f 
8821 
eut 

9441 


i  ttmmmmmt  Diagnostics  tmmtmtmtmmmtmmm* 
i  muuttinuatmmtiimmuuummumutmttumttttnmm 


i«propep  oxides  present 


9661 
8679 
B680 
8679 
8701 
8710 
8721 
87  .SO 
8740 
8750 
8760 
8770 
8780 
8790 
8800 
8810 
8820 
8830 

am 

8851 

8860 

8870 


Alert2 i  ! 

PRINT  "Error  in  OXIDES" 

COSUB  Son 
STOP 

OT0  Tens 
i 

PRINTER  IS  7,1! 

GOSUB  Spec 
COSUB  Prtflg 

PRINT  LIN(1)/T  =  *,T,LIN(1) 

1=1 

pRINT  USING  Inign;FcKI),Nd).Fe»d+l).Nd+i>.F»i(l+2)  NU+2).Fi6(I+3),NTl*3) 


Sun: 


If  ac'nated  b?  lil,  orints  after  each  itergtien 


PRINT  USING  IjwuniFeKD.N'D.Fotd+D.NfWl/eKWl  NfWl.FetfWI.NfW) 

i=10 


PRINT  USING  l«5«n;Fg6d),Nd),Fo6d*i),Nd+l),Fo6iI+2),N(l+2).Fei(I+3) , NT  1+3) 


PRINT  USING  IfwgAiFiBdi.Ndl.Frid+O.NU+D.FoiU+Ei.Nd+ai.FeFlMJ.NU+S) 


PRINT  USING  tMSun:.Fo*d>,Nd>,T!i*UH>,N(I+i>.Fi»d+2'  80+2)^06(1*3) ,N'I+3) 
I*?? 


PRINT  USING  In*il«iFl4(I)fNd)lFo4(I+l),N(I+l),Ft4(l+2)  N(IJ2),Fl4(I+3),N(H3) 

4"" 


PRINT  USING  Insg«)Fo»dt,Nd).Fg4'I+l),N(!+l) 

PRINT  USING  Insun; *TiC\ tic, "fret*  Net, lTiO‘ ,Til,*T.233", 
PRINT  USING  Ins».ii,Ti30^,n'!,"Tid2’,fj4  "TiH'.fin  "Cr"> 

lUknr  lift*  u-r  hnlr  v  •  1 


8891 

8900 

8910 

8920 

9931 

8940 

0950 

8960 

0970 

89811 

8991 

9000 

9010 

9020 

9031 

9040 

9051 

9060 

9071 

vieo 

9090 

9lDD 

9110 

9120 

9130 

9140 

9150 

9160 

9170 

91B0 


.  .  _n  ■  WWArrv  4niV.li  MAU  l 

Insun:  IHAGE  J(8A,IIZ.i>ME; 3X> 

Sunl:  Slm=FNSti+Tin+Tic+Net+FNTsuni 

PRINT  LINUL'Ti  BALANCE  ".Nti.Sun 
3iin=FNS«+Tii+38Ti2+54Ti3*2*Ti4 
PRINT  “0  BALANCE  VNea.Sgn 
S«n=Tin+26N2+Ne+Hcn+Hnc«+Cn*C2n+Cnn 
PRINT  "N  BALANCE  "  Nn,S«n 
Sin=Nc--Ncc+C»+Ci2+l lc+Gr 
PRINT  "C  balance", Nc , Suit 
S'jn=Hh*Y+2«Dhir*y 
PRINT  “H  balance", Nh .Sun 
°RINT  USING  'KfK‘;*Acc=  ",Acc 


Tl.I 

Sr 


Aopeo+s  in  find  autcut 


PRINT  USING  "K.K »//".'  "Atsim  =  ",Atsun 

’PINT  LINI2) ; T«Bv£0 > , "Test  for  tpvlliUrivn'.LINd) 

hr. hh.,,.  «...  ........  ...  ,  if  •,TAI.(7»1.' 


PRINT  SPAd),  K*;TAB(10);"IT«nef  ic";TAB(30)j 
PRINT  "Fijnctiotial((=Nuneric)".lIN(ii 
FLOAT  3 

IP  N9T  loiflaq  AND  NOT  Lflaq  THEN  A  +0=At l=ftt2=At3=A*4=l 
PRIM  "Ti<g>";TABll(!),Atil/KhjTAB(30Vi 
PRINT  *Ti0';TAB(H),An/Xl!TA?(30!,Ti»X 


PRINT  ”Ti203*,TABllO).,SQR(At2)/K2 '.TABI31I ) ,Ti  extSBRd i 
PRINT  •Ti3O5",TABtl0).SOR(At3)/X3;TAB(30),Ti8XtX*S8R(Ti»X) 
PRINT  Ti02  jTABilO' .At4/K4:TAB(3l),Ti802 


PRINT  ’TiN*;TAB(10),i/Ktin:fAB(30),fi<2 


PRINT  "nC'iTABflOM/KticjTABTJOr.TiSAcc 
IF  Liuflag  OR  Lflaq  THEN  PRINT  LINT! ) 


IT  Lixflag  08  Lflaq  JHEN  PRINT  USING  Issuh; ‘gTNet >*.ftrS. "adiO! " .At* ,"a(Ti203)" .fit? 
IF  Lcxflaq  OR  Lflaq  THEN  PRINT  USING  iNsuni'addQS)'  Atj/adiOE)",^ 

L.  Actmtm  in  liquid  nide 

STAnMRO 

PRIkTE*  IS  7,1 
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9200 


I  HUUUUttltMUIIMIllMItilllllllUtOUtllttllllintMIIIMtmttltlM 

1 


9238 

9248 

92S0 

9268 

9278 

9281 

9290 

9300 

9311 

9320 

9330 

9340 

mi 

9368 

»37l 

9380 

9391 

9400 

9410 

9420 

9431 

9448 

9458 

9460 

9471 

y4BII 

9490 

9510 

9511 
9520 
9538 
9541 
«550 
9568 
9571! 
9580 
9598 
9600 
9610 
9620 
9630 
')640 
9650 
9660 
9670 
9680 
9691 
9700 
9710 
9720 
9738 
9740 
9750 
9760 
9770 
9780 
9790 
9800 
9810 
9820 
9830 
9840 


i  1 


Ti 

no 

T102 

Tl+ 


At 

CO 

C02 

H 

OH 

% 


(nit  uttd] 


H20 


init  mid) 
Inot  uid 


& 


0 
02 

T1CO) 

Tins) 

Till) 

Tm' 


CN 

ON 

C2H 

C2N 


51 


HCN 

HNCO 


Inal  usidl 


C2H2 

C3 

TiO(l) 

TiO'f) 

Ti20J(l) 

T i203 1  s' 

Ti305!l) 

U305IS' 

Ti02U) 

U02(s' 

riNd) 

TiK(s) 


lie 

Tit 

Hit 

Hit 

Or 


'll 

hi 

TiC 

ll? 

in 

li4 

7i4 


■  in 
Tin 


flout 

Flaql) 

CflSa 


ri$ 

FlagZ 

Fhg2 


F  aq3 
Flag3 


lag3 

Flog4 

Fluo4 

Tinflaq 

Tinflaq 


tmtuttt 


Ktir 
K  tie 
K 1 1 
till 


SI 

*1 

V2=?BK(Knj 
K?>SWO'n  > 
H>SBIKk'n' 
KMMtKN' 
K4 
M 
K'in 
K  tin 


Loiflau  ir  Iflaq  indicates  pnsnnce  of  cent lr>itois  liquid  ohoti 
A10,  At i,  At2,  At3,  At4  far  activity  Udoii  silytiin' 
of  Hit.  ^0,  fi203,  T1305,  T1C2 
Act  =  artmty  if  carbon 
X  =  509(02) 

Y  =  S»(H2) 

2  =  S8R(H2> 


Etlt 


Oanafcles  usto  is  Niytsn's  nttnid  colcilatii'1: 

No  mtql  or  oiidi  X,Ti 

One  nidi  inly/nital  inly  t  <Ti  iF  T|N  *  1  sal  if  otldi,  no  natal ) 
Two  audis/Hital  *  1  oxide  - 

lp  oil  casts,  Z  and  Acc  are  the  last  two  oaniahUs  when  gpprooriate 
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